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In this review paper, we discuss the influence of weak quenched disorder on the criti- 
cal behavior in condensed matter and give a brief review of available experimental and 
theoretical results as well as results of MC simulations of these phenomena. We con- 
centrate on three cases: (i) uncorrelated random-site disorder, (ii) long-range-correlated 
random-site disorder, and (iii) random anisotropy. 

Today, the standard analytical description of critical behavior is given by renor- 
malization group results refined by resummation of the perturbation theory series. The 
convergence properties of the series are unknown for most disordered models. The main 
object of this review is to discuss the peculiarities of the application of resummation 
techniques to perturbation theory series of disordered models. 



1. Introduction 

There are two main questions we want to discuss in this paper. The first question is 
a physical one: how does disorder influence criticality? Taken that a non-disordered 
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( "ideal" ) system possesses a phase transition and a critical behavior governed by 
some scaling laws, will these universal laws be altered by disorder introduced into 
the system? The second question elaborated below is related to the unknown con- 
vergence properties of the series that result from the perturbative treatment of the 
corresponding models. 

The first question we address here is rather specific since it concerns properties 
in a narrow region of the phase diagram in the vicinity of a critical point. Nonethe- 
less, in order to give an extensive description of current studies in the field we are 
forced to restrict the subject further. In particular, in the present paper we will be 
interested only in static critical properties leaving the bulk of dynamic phenomena 
outside of our scope. We will choose three-dimensional systems of infinite size ignor- 
ing the infLuence of surfaces on the critical behaviour though such problems are very 
attractive 1^1. Furthermore, we will address the subtle limit of weak disorder where 
the random system can be treated as perturbed with respect to an "ideal" one; this 
disregards the class of phenomena relevant in the strong disorder limit, like those 
appearing in the vicinity of a percolation threshold. Finally, we will touch only the 
quenched "frozen" disorder. To summarize, the main scope of our review will cover 
the static criticality of random quenched three-dimensional systems. In spite of our 
restriction of the subject, a lot of fascinating and still unsolved problems remain as 
we will try to convince the reader in our subsequent exposition. 

To get a quantitative description of the phenomena occurring in the vicinity of 
the critical point it is standard to rely on renormalization group methods 0. In this 
approach the physical quantities describing the critical behavior are obtained in the 
form of perturbation theory series. For some ideal (undiluted) systems these series 
have proven to be asymptotic and to have a zero radius of convergence. Therefore, 
special resummation methods are used to extract reliable data from the series. How- 
ever, so far there is no proof of the asymptotic nature of the perturbation theory 
series for disordered (diluted) models. Moreover, there exist serious indications of 
possible non-asymptotic divergences. Nonetheless, the bulk of results that were ob- 
tained from the perturbation theory expansions of disordered models are presented 
as if they were asymptotic. Here arises the second question we address in this paper: 
how does the possible non- asymptotic nature show up in the analysis of perturba- 
tion theory expansions for the models with weak quenched disorder? Putting the 
first and the second questions together, the main goal of the current review will be 
to discuss what information about the influence of weak quenched disorder on the 
critical behavior can one extract working with renormalization group perturbation 
theory series. 

The set-up of the paper is as follows. In the next Section ^ we discuss different 
models of weak quenched disorder and introduce the physical quantities used for 
their description. Here, of main interest for us will be systems with non-correlated 
point-like defects, systems with long-range correlated disorder and systems with 
disorder in the form of random anisotropy. In Section we give examples of ex- 
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perimental realizations of such systems. These include: site-diluted magnets, He"* 
in porous medium, magnets with extended impurities, polymers in random media, 
some amorphous rare-earth-transition metal alloys. The effective Hamiltonians for 
the disordered models under consideration are given in Section ^ together with the 
main relations of the field theoretical renormalization group theory. In Section ^ we 
describe the methods of resummation of divergent series and in Section ^ we dis- 
play the results obtained by their application to the field theoretical renormalization 
group expansions. We conclude and give an outlook in the Section m. 



2. Different types of weak quenched disorder 

We start our consideration from one of the simplest and most carefully studied 
models which is introduced as a generalization of the hierarchy of the lattice Ising 
model, the XY-model and the Heisenberg model. Consider a simple (hyper)cubic 
lattice of the dimension d, to each site of which it is prescribed a m-component 
vector S with a fixed length (for convenience of further notation one usually sets 
l^l = ^/m). Imposing a pair interaction with the energy proportional to the scalar 
products between pairs of spins, one defines the Stanley model (it is also known as 
m-vector, 0(m)-symmetric or generalized classical Heisenberg model). The Hamil- 
tonian of the Stanley model reads i: 

R,R' 

where the d-dimensional vectors R span the sites of a simple (hyper)cubic lattice 
(see Fig. 0) and J(|R — R'|) is an isotropic and translationally invariant short-range 
interaction. The Stanley model has further been generalized by analytic continua- 
tion in TO to describe critical properties of self-avoiding walks (to — 0), spherical 
ferromagnet (to = oo), the Gaussian model (to = —2), etc. 

A many-body system described by the model (^ can undergo a continuous 
phase transition into a ferromagnetic phase. Yet, the transition occurs only for lat- 
tice dimensions d larger than a certain lower critical dimension dLow The numerical 
value of c/low is a function of the spin component number to. If the Hamiltonian (|^) 
remains invariant under continuous rotations in the space of vectors S, i.e. to > 2, 
one can show that dLow = 2 Otherwise, for to = 1 (the Ising model) the continu- 
ous 0(to) symmetry degenerates into a discrete one corresponding to the invariance 
of the Hamiltonian^(^ under discrete fiip of all spins to the opposite direction. In 
this case dLow = 10. 

In the theory of critical phenomena one is interested in the peculiarities of the 
behavior in the vicinity of the critical temperature Tc > 0. The critical temperature 
essentially depends on microscopic details of the model. However, asymptotically 
close to Tc the thermodynamic and correlation functions for very different systems 
may depend on the external parameters by the same power-like functions. These 
functions are primarily characterized by their critical exponents. 
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Fig. 1. The Stanley model describes a system of m-component vector "spins" occupying sites of 
a d— dimensional (hyper)cubic lattice and interacting by means of a short-range translationally- 
invariant force. 



The critical exponent a; of a physical observable 0{t) asymptotically close to 
the critical point Tc is defined by i: a; = lim^^o ^"ifj^ i where t — (T — T^)/T(. is 
the reduced distance to the critical point. For instance, the magnetic susceptibility 
X diverges as: 

X~r±|r|-^ r->0 (2) 

where 7 is the susceptibility critical exponent, while F-i- and r_ are critical ampli- 
tudes above and below the critical point respectively. The power law of type (||) is 
exact in the asymptotic regime r ^ 0. 

The definitions of the most common static critical exponents of a magnet are 
provided in Table |l|. 



Table 1. Critical exponents governing the power-like behavior of thermodynamic and correlation 
functions asymptotically close to Tc- Here, M is the magnetization, x is the magnetic susceptibility, 
Ch is the heat capacity at constant external magnetic field H\ {5M{fi)6M{K)) is the pair correlation 
function of the magnetization fluctuations and J is the correlation length. 

M (X (-r)'3(T ^ T", H = Q) (5M(0)5M(R) > oc exp(-|R|/0 (T ^Tc, H = 0) 

H oc |M|*sgn(Af) (T = Tc, H 0) ^ oc t-" (T T+ , H = 0) 

X oc T-'' (T ^ T+, H = 0) (<5Af(0)<5A/(R) > oc |R|-(''-2+'') (T = Tc,H = 0) 

Ch oc r-° (T ^ T+, H = 0) 



In contrast to the critical temperature, the values of the critical exponents do not 
depend on microscopic details of the Hamiltonian and are determined only by the 
global properties of the model. For a short-range interaction the global properties 
are the lattice dimension d, the dimension m and other symmetry properties of 
the order parameter. Due to the fact that the critical exponents may be identical 
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for systems with very different microscopic nature, the critical exponents are called 
universal. Subsequently, if the critical properties of two systems are described by 
the same set of scaling functions and an identical set of critical exponents they are 
said to belong to the same universality class. 

In the case of the 0(m)-symmetric model (Q) the critical exponents are functions 
of the lattice and spin dimensions d and m. However, if d is greater than some upper 
critical dimension dup = 4, the critical exponents actually do not depend on m and 
take mean field values 7 = 1,1^= l/2:'7 = 0, /? = 1/2. For dLow < < ^Up the 
critical exponents show non-trivial dependence on the spin dimension m and need 
further consideration. In particular, for space dimension d = 3 the critical exponents 
of the model (|^) have been calculated so far with very high accuracy. The standard 
values of the critical exponents are listed in the Table M. 



Table 2. The stanrfprd numerical values of the critical exponents of the d = 3 Stanley model 
(results of the Ref. 0). The cases m = 1, 2, 3 correspond to Ising, XY- and Heisenberg models, 
m = reconstitutes the so-called polymer limit. 



m 


7 


V 


»? 


/3 


Q 





1.1596 ±0.0020 


0.5882 ±0.0011 


0.0284 ±0.025 


0.3024 ± 0.0008 


0.235 ± 0.003 


1 


1.2396 ±0.0013 


0.6304 ± 0.0013 


0.0335 ±0.025 


0.3258 ± 0.0014 


0.109 ± 0.004 


2 


1.3169 ±0.0020 


0.6703 ±0.0015 


0.0354 ± 0.025 


0.3470 ± 0.0016 


-0.011 ± 0.004 


3 


1.3895 ±0.0050 


0.7073 ±0.0035 


0.0355 ± 0.025 


0.3662 ± 0.0025 


-0.122 ±0.010 



The model (|l|) corresponds at criticality to a wide range of real systems, not 
always magnetic, and it is a convenient starting point to study critical phenom- 
ena. However, it suffers from obvious restrictions. Typical for real substances are 
properties which can be considered as deviations from the "ideal" structure of the 
Hamiltonian (j^). Various types of defects spoil the strict translational invariance 
of the inter-spin interaction J and/or can be effectively represented by additional 
terms in the Hamiltonian. In both cases, disorder introduces stochastic (random) 
variables in the Hamiltonian of the system. In the following, we will treat random- 
ness introduced by either site dilution or random anisotropy. 

The presence of non-magnetic impurities can be modeled by a class of site- 
diluted models Q. A site-diluted Stanley model is introduced by the Hamiltonian: 

where cr is an occupation number: it equals 1 when the site is occupied by a 
magnetic atom and otherwise (see Fig. ||). Geometrically, the vacancies — ^ 
can be distributed independently according to the probability density 

P(cj^) - (l-p)^(cfi)+p<5(l-cfl) (4) 

with p being the concentration of magnetic atoms. In the absence of correlations 
between the vacancies the ordered phase of the model (||) persists with a decrease 
of concentration p up to some value Pc- The latter depends on the lattice type and 
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Fig. 2. The weakly diluted Stanley model is characterized by a small amount of randomly 
distributed (quenched) non— magnetic atoms which spoil the translationally-invariant structure of 
the spins in Fig. 



dimension and is equal to the percolation threshold ppac a. Here, we are not inter- 
ested in phenomena occurring in the vicinity of the percolation threshold p ~ Pporc 
but rather we address the limit of weak dilution p ~ 1. Since in the critical region 
impurities cause, in particular, a shift of the critical temperature, this property 
can be taken as a genuine property of disorder due to non-magnetic defects. In 
turn, the scale of the correlation of non-magnetic impurities becomes apparent by 
the fluctuations of the local temperature of the magnetic phase transition. In what 
follows below we will not consider frustrations in the exchange interaction J{R) in 
which may lead to the spin-glass phase as far as the main problem of interest 
for us is the second order ferromagnetic phase transition. 

From the viewpoint of dynamics one may distinguish two main types of disorder. 
If the dynamics of the impurities has a characteristic time that is comparable to the 
relaxation times in the pure system, the impurity variables are treated identically to 
the "pure" dynamical variables and are a part of the disordered system phase space. 
The corresponding annealed disorder □ has been subject to numerous investigations 
and reveals trivial results at criticality. The so called Fisher renormalization 113 
states that when a heat capacity critical exponent of an undiluted system a^ure is 
positive, then the critical exponent a of the corresponding annealed diluted system 
changes to a = and any other critical exponent x of the annealed system 

is determined by the corresponding pure one {xpure) by a simple normalization 
of the form: x — /y^''' . This explains the prevailing interest in systems with 
quenched disorder where the impurities can be considered as fixed. To obtain the 
thermodynamics of the quenched system one needs to perform the configurational 
average over an ensemble of disordered systems with different realization of the 
disorder. 

The influence of uncorrelated weak disorder of type (jj) on the critical behavior 
of the model (||) can be predicted by the so-called Harris criterion til. It states that 
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only if the heat capacity critical exponent Ofpuro of the pure (undiluted) system 
is positive, i.e. the heat capacity diverges at the critical point, quenched disorder 
causes changes of the critical exponents. From Table I one concludes that only the 

tng model is affected by short-range-correlated quenched disorder at criticality 
Moreover, the heat-capacity critical exponent of the weakly-diluted Ising model 
should be negative as is suggested by the inequality v > 2/d proven for diluted 
systems 0. 

To explain possible generalizations of the model with site occupation dis- 
tribution let us note, that (as it will be seen from Section ||) for well behaved 
distributions the two first moments as far as they exist determine the critical be- 
haviour of the model. For the distribution (^ one gets: 

(cr) = p, (5) 

g(|R- R'l) = (CRCR,) - (CR>2 = p(l -p)5r_r,, (6) 

where (. . .) means averaging with the distribution function (^), and (5r_r' is Kro- 
necker's delta. 

For the pair correlation function (^ one may consider Il3 more general cases. 
In the first of the two generalizations treated in this review the quenched random 
impurities are strongly correlated in Ed dimensions and randomly distributed over 
the remaining d~ Sd dimensions while in a second generalization we will treat long 
range correlations as described below. In the first case, the impurities are envisaged 
as parallel e^-dimensional objects, each extending along the coordinate directions 
symbolized as R| | , whereas in remaining d — Sd dimensions symbolized as R_l they 
are randomly distributed. The pair correlation function in this case is defined to 
yield 0: 

5(|R-R'|) = i'o'^''~'"(R± withuo>0. (7) 



The case of point-like uncorrelated impurities (and the pair correlation function 
is reconstituted from (|^) by — 0. Systems with impurities distribution governed 
by the pair correlation function (^ are not covered by the aforementioned Harris 
criterion. To show this, let us note that due to fluctuations in concentration appears 
the local transition temperature will vary. Denoting the deviation from the local 
transition temperature by t{R) one can determine the mean deviation r^d to the 
transition temperature in a region of correlated spins, i.e. in a volume with 
correlation length ^, as the average of t(R) over this region. 



1 / d'^i?r(R). (8) 

4 Jf'' 



For convenience, we pass to the continuum space R with subsequent change of 
Kronecker's deltas to (5-functions and of sums to integrals. 

The mean values (|^) vary from one correlated volume ^ to another. One can 
only speak consistently of a phase transition in the whole sample at a precise tran- 
sition temperature if these variations become negligibly small as the transition is 
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approached. The variation of the mean values rjd in different correlated regions is 
characterized by their variance 

= i / ^'rI d'R'gi\R-Il'\). (9) 



For the weak impurity concentrations the shift in the local transition temperature is 
proportional to the concentration. This allows one to substitute the pair correlation 
function in by the impurity pair correlation function (^. Evaluating (||) and 
using ^ ^ r^^p (with Vp being the correlation length critical exponent of the pure 
model) one finds that for t ^ 0: 

This tends to zero, i.e. the variance becomes negligible, for t ^ if the exponent 
of T is positive; if it is negative, however, the variance A eoes to zero more slowly 
than T itself, and the mean value r becomes meaningless 113. This is the indication, 
that the defects qualitatively change the transition, and it leads to the generalized 
Harris criterion for the relevance of extended defects on the critical behaviour: 

ed>d~-. (11) 
i/p 

Putting in ( [Tl| ) £dj-^ one recovers the Harris criterion for the point-like weak 
quenched disorder til. 

The second generalized correlation function that^e treat here describes a power- 
law decrease of the correlations at large distances O: 

5(|R-R'|)^|R-Rr"- (12) 



In what follows, defects which are correlated according to (|12|) with a > d will 
be called short-range correlated, and correspondingly, with a < d - long-range- 
correlated. The correlation function ( p^ ) has a direct physical interpretation: straight 
lines of impurities of random orientation are described by the case^a = d— 1, whereas 
randomly oriented planes of impurities correspond to a — d~2i3. Moreover, differ- 
ent non-integer values of a may be interpreted as randomly distributed impurities 
of fractal nature O. Weinrib and Halperin ta showed that in the presence of long- 
range power-law correlations in the disorder the Harris criterion is also modified: 
for a < d the disorder is relevant, if the correlation length critical exponent of the 
pure system obeys u < 2/a, while for a > d the usual Harris criterion for the in- 
fluence of short-range-correlated point defects applies. The derivation is similar to 
that presented above for the case of extended impurities. 

As mentioned above, apart from site dilution as a second example of randomness 
we treat a random anisotropy present at all sites. Let us introduce a model for 
Quenched weak random anisotropy proposed by Harris, Plischke, and Zuckermann 
lla. This model describes a regular lattice of magnetic ions, each of them being 
subject to a local anisotropy field of random orientation in the form of an easy 
axis. The magnitude of this field is the same for all sites. The Hamiltonian of the 
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random anisotropy model (RAM) then has the form of the Hamiltonian (|l|) with 
an additional anisotropy term: 

n = -^Yl Jm-^'\)SR-Sti'-DY,i^-[iSnf. (13) 

Here, > is an anisotropy strength (the case of _D < corresponds to the 
presence of easy planes in the model), and zr is the unit vector pointing in the 
local (quenched) random direction of the uniaxial anisotropy. 

Strong disorder at large ratio D/J, should destroy the second order phase tran- 
sition. For small D/J this is not so obvious. In addition to the global variables of 
a regular magnet (i.e. lattice dimension, type of interaction and spin symmetry) 
the low-temperature ordering in RAM is influenced also by the distribution of the 
random variables 2 = zr, in (|l^). For non-correlated z the low-temperature order- 
ing depends on the probability distribution p{z) of the direction of anisotropy on 
a single site. As introduced by Aharony 111 two kinds of distribution functions that 
are relatively simple to study and that correspond to real physical situations, are 
generally considered. The first one corresponds to the isotropic case, when the ran- 
dom vector z points with equal probability in any direction of the m-dimensional 
hyperspace 

with T(x) being Euler's gamma-function. The second one corresponds to the so- 
called cubic anisotropy, when the vector z points to any of the 2m directions ki 
along the edges of the m-dimensional hypercube: 

P(^) = ^ Y.^S^'^'^z - h) + <5(")(f + k)]. (15) 
2m ^ — ' 

1=1 

The rationale for such a choice is to mimic the situation when an amorphous magnet 
still "remembers" the initial (cubic) lattice structure. Other distributions may be 
considered as well. 

For the isotropic case (|l^) arguments similar to those of Imry and MaE^ applied 
to the random-field Ising model demonstrate the absence of magnetization below 
d = A. Applying domain arguments EiJ one supposes that the ground state will break 
into finite clusters of linear size L. This is possible if the bulk energy DL~'^/'^ 
gained is larger than the domain wall energy. For isotropic spin systems the latter 
energy is ~ JL^^^. This leads to the conclusion that below four dimensions a break 
up into domains is favored with a characteristic domain size Lq ~ {J/D)'^/^'^^'^\ 
Other arguments were applied by Pelcovits et al. El. Assuming a ferromagnetic state 
they calculated the correlation function of transverse magnetization and obtained 
its divergence for < 4. However these arguments do not consider the anisotropic 
distributions p{z). 
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In the next section we introduce various physical systems described by disor- 
dered quenched models (||) , (^3|) . A brief review of available experimental and MC 
data for these systems is given as well. 

3. Realizations of weak quenched disorder: experiments and 
Monte Carlo simulations 

In this section, we review critical phenomena in substances and their computer 
simulation that can be considered as realizations of the models discussed above. We 
start from weakly site-diluted quenched Stanley models. Among these, the bulk of 
results is known for the weakly diluted quenched Ising model as is suggested by the 
Harris criterion. 

3.1. Site-diluted uniaxial antiferromagnets 

Typical experimental realizations of the weakly diluted quenched Ising model are 
given by crystalline mixtures of two compounds. The first is an "Ising-like" ani- 
sotropic uniaxial antiferromagnet with dominating short-range interaction (e.g. 
FeF2, MnF2), the second one is non-magnetic (ZnF2). Mixed crystals (FepZni_pF2, 
MnpZni_pF2) can be grown with high crystalline quality and very small concentra- 
tion gradients providing an excellent realization of random substitutional disorder 
of magnetic ions (Fe"*"^, Mn+^) by non-magnetic ones (Zn+^). 

Already the first experimental studies gave evidence of new critical behavior 
caused by a weak quenched dilution. In nuclear magnetic resonance measurements 
of the magnetization in Mno.864Zno.i36F2 the value of the magnetization expo- 
nent (3 = 0.349 ± 0.008 was found to differ strongly from that in an undiluted 
sample. Within a few years this result was corroborated by nuclear scattering mea- 
surements of the_magnetic susceptibilitv_and correlation length critical exponents 
in FepZni_pF2 cilta and MnpZni_pF2 c3 at different dilutions. In particular, the 

Sbter study showed that critical behaviour of the susceptibility in Mno.75Zno.25F2 
is governed by a critical exponent 7 which, consistently with the Harris criterion, 
differs from that of the "ideal" magnet (see Fig. |3|). The linear birefringence mea- 
surements revealed the cusp-like behavior of the specific heat at the transition point 
with an exponent a = — 0.09± 0.03 EJ. This experimentally proved that within the 
error bars the hyperscaling relation dv + a = 2 is satisfied. 

Crossover phenomena in diluted systems are governed in addition to the tem- 
perature by the concentration of the magnetic component. The experimentally ob- 
tained exponents are often reported to be effective ones (i.e. temperature- and 
dilution-dependent). However, already in the first experiments on the critical be- 
havior of quenched diluted Ising-likc antiferromagnets it appeared possible to reach 
the asymptotic region. Thus, studying the critical regime in Ref. Efl the authors 
found neither a region in reduced temperature r where one finds "pure" Ising ex- 
ponents nor any evidence of crossover from pure to random exponents. This was 
explained by the crossover either taking place outside the critical region or being 




Fig. 3. Neutron scattering measurements of the susceptibility x(0) in Mno.75Zno.25F2 Ej. The 
solid lines are the results of fits to single power laws with exponents 7 ~ 1.364 above and below 
Neel temperature Tjy. The critical behaviour is governed by a simple power law in the reduced 
temperature interval 4 ■ 10"'' < |r| < 2 ■ 10~'. 



too slow. In Ref. the crossover from pure to diluted critical behavior was studied 
and a magnetization exponent P — 0.36 ± 0.01 was found that does not change 
under dilution for p < 0.05. The cross ove r occurs within a very narrow range of r 
at relatively large values of r. In Ref. c3 excellent agreement of the measured ex- 
ponents 7 = 1.364 ± 0.076 and ly — 0.715 ± 0.035 with the theoretical values of the 
diluted Ising model was obtained for the temperature range 4-10"'* < r < IQ-^ and 
p — 0.75. The value of the critical exponent (3 was also the subject of a crossover 
analysis in Refs. Ref. El concludes that the experimental errors are too large 

in order to distinguish between the pure Ising model and the diluted Ising model 
critical behavior. In Ref. El no crossover was found after the correction-to-scaling 
had been taken into account, and the diluted Ising model critical behavior was 
found in the whole temperature range. 

It is known that the diluted Ising magnet in a uniform magnetic field H along the 

Biaxial direction exhibits static critical behavior of the random-field Ising model 
. Such experiments give additional information about the critical behavior of the 
diluted Ising model when performed for H ^ Oc3. A comprehensive account of 
the experimentally measured valuea_pf the critical exponents of diluted uniaxial 
magnets is given in a recent review t3. 
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3.2. Disordered Heisenberg-like magnets 

Another experimental realization of the weakly diluted Stanley model are given 
by magnetic systems described by the site-disordered Heisenberg model Some 
amorphous magnets as Fego+xZrio-x, Fega^MyZrio ( M = Co, Mn, Ni) E3o and 
transition-metal basecLmagnetic glasses EaE3 as well as disordered crystalline ma- 
terials as Feioo-xPtx FeyoNiso E3 and Eu-chalcogenide sohd solution systems 
(EuxSri_xSo.5Seo.5, EuxSri^xS, EuxBai_xS, EuxLai_xS) ^"''^'^ belong to this group 
of systems. 

According to the Harris criterion the critical behaviour of the diluted Heisenberg 
model should remain unaltered in its asymptotics with critical exponents of the 
c? = 3 pure Heisenberg model (see the last line of Table However, the experiments 
give a wide scattering of the values of the critical exponents. Often the reason is 
that the asymptotic region still is not reached and only the effective exponent is 
observed. In particular, the value of the magnetic susceptibility effective critical 
exponent 7e/ / with decrease of distance to the critical point goes through a peak to 
its asymptotic value which coincides with magnetic susceptibility critical exponent 
of pure Heisenberg magnets. The dependence of 7e// on reduced temperature r is 
presented in Fig.|| taken from the recent paper a. The effective behaviour also was 
studied by Monte Carlo simulations lj, the results show that 7e// has a maximum 
while approaching the asymptotic regime. Such form of temperature dependence 
of effective su^eptibility critical exponent is corroborated by recent theoretical 
investigations c3. 

1.8 
1.6 

1.4 

1.2 

0.01 -J 0.1 

Fig. 4. Experimentally measured isothp»nal susceptibility effective critical exponent "feff for 
Fego-xMnxZno for different values of xM {t = (T - Tc)/Tc). 




3.3. He* in a porous medium 

The listed experimental evidences of the critical behavior of the weakly diluted 
quenched Ising model confirm the Harris criterion. However, the bulk of the data 
was obtained for the magnetic phase transition. Due to the universal character of 
the critical exponents the model (El) can describe universal properties of critical 
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phenomena of completely different microscopic nature. An example is the transi- 
tion of liquid helium-4 to the superfluid state which belongs to the Stanley model 
universality class at m = 2. The uncorrelated point-like disorder (Q) in the case 
of a fluid can be introduced by a porous medium cil. Recently, high precision ex- 
periments on the critical behavior of liquid helium-4 near the superfluid transition 
in a porous medium e3 confirmed the irrelevance of quenched disorder, again in 
correspondence with the Harris criterion|_since the specific heat exponent near the 
superfluid transition is slightly negative e3. More data|_about effect of disorder on 
the A-transition of He^ may be found in recent review E3. 



3.4. MC studies of random Ising model 

Another way to study the critical behavior of the model (||) is to use computer 
Monte-Carlo experiments. However, in contrast to experiments, the first studies 
of critical behaviour of the weakly diluted Ising model on a simple cubic lattice 
l3 revealed critical exponents that were identical within the numerical error with 
the corresponding exponents of the pure (undiluted) system within a wide dilution 
region. Later, these data were objected by subsequent MC simulations on larger 
lattices c3, that obtained critical exponents varying continuously with the mag- 
netic sites concentration p. Indications of a change of the order parameter critical 
exponent (3 upon dilution initiated an extension of the studies to determine the 
other critical exponents and to check the scaling for disordered systems. Finally, 
the Swedsen-Wang algorithm that avoids the critical slowing down of the relaxation 
was applied to the d = 3 diluted Ising model Ezl and resulted in critical exponents 
for the susceptibility and correlation length that were independent of concentration 
over a wide range_DLdilution. i— ii— i 

Due to Refs. Ellea and especially EaE3 it became clear that the concentration 
dependent critical exponents found in previous MC simulations are effective ones, 
characterizing the approach to the asymptotic region. The effective exponents 7, /? 
and C = 1 — P (the last one describes the divergence of the magnetization-energy 
correlation function) were shown ca to be concentration dependent in ihe concen- 
tration region 0.5 < p < 1. These data were refined three years later e3 resulting 
in more accurate estimates for the above mentioned exponents and the critical 
exponent i> of the correlation length with continuously varying values. 

The iritical behavior of the d = 3 diluted Ising model was reexamined recently 
in Ref. O. The study was based on the crucial observation that it is important 
to take into account the leading correction-to-scaling term in the infinite volume 
extrapolation of the MC data. The simulations confirmed the universality of the 
critical exponents of the d — 3 diluted Ising model over a wide region of concentra- 
tions (see Fig. |^) . In particular the value of the correction-to-scaling exponent uj 
was found to be w = 0.37 ± 0.06 which is almost half as large as the corresponding 
value in the pure d — 3 Ising model to — 0.799±0.011 i. The smallness of lo in the di- 
lute case explains its importance for an analysis of the asymptotic critical behavior. 
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The other exponents were determined as: /3 = 0.3546 ± 0.0028, 7 =p-l|442 ± 0.010, 
V = 0.6837 ± 0.0053. We refer the interested reader to the papers E3e3 where the 
numerical values of the exponents obtained in the MC simulations are listed. 




0.0 0,1 0,2 0.3 0,4 0,5 



Fig. 5. Determination of the correction-to-scaling exponent uj in the MC simulations of Ref. a. 
Quantity Qd^^ = 2^+^/'^ is plotted for different dilutions 0.4 < p < 1 and lattice sizes 8 < L < 128. 
One can see that at p = 0.9 the system is crossing over from the pure Ising fixed point to the 
diluted one even for L = 128. The solid lines correspond to a fit a; = 0.37 yielding the same infinite 
volume extrapolation for all p < 0.8. 



3.5. Magnets with long-range- correlated quenched disorder 

Whereas at present many experimental and MC simulation data are available for 
magnetic systems with random-site uncorrelated disorder, the critical behavior of 
systems with long-range-correlated disorder is much less investigated. Therefore, 
in this subsection that is devoted to systems with this type of disorder we present 
only a few works where its influence has been studied. 

In magnets, long-range-correlated disorder may be present in the form of con- 
tinuously distributed dislocations and disclinations, so-called extended structural 
defects. These defects may have the form of lines or planes of random orientation E3 
or may form some sponge-like fractal objects, which are considered as aggregation 
clusters 113. The particular case of systems with dislocation lines or planes of parallel 
orientation, is investigated in Ref. P^, 

Recently, Ballesteros and Parisi Ej presented Monte-Carlo simulations of the 
diluted Ising model in three dimensions with extended defects in the form of lines of 
parallel orientation. Finite-size scaling techniques were used to compute the critical 
exponents of these systems, taking into account the strong scaling corrections. The 
obtained value for the exponent v is compatible with the analytic predictions 0. In 
Ref. the three-dimensional so-called thermally diluted Ising systems were studied 
in Monte-Carlo simulations. Here, at the characteristic ordering temperature 9, the 
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spins of the dominant type (concentration p > 0.5) are taken as the location of the 
magnetic atoms, and the rest are taken as magnetic vacancies. The structure of the 
realization is fixed thereafter for all temperatures at which the magnetic interactions 
are subsequently investigated. If happens to coincide with the critical temperature 
of the undiluted Ising model, the vacancies are in randomly located points, but 
with a long-range-correlated distribution. It was found that the simulated critical 
exponents of thermally diluted Ising systems at criticality agree fairly well with 
theoretical predictions of Weinrib and Halperin 0. 

We adduce here also the interesting result of Monte-Carlo study of three-dimen- 
sional Ising model with bond disorder, which is random in the vertical direction 
and correlated in = 2 horizontal directions The results demonstrate a sharp 
second-order phase transition, with exponents of new universality class. However, 
it seems not to be of the same universality class as the model of extended defects 
presented by Dorogovtsev Il3. 



3.6. Polymers in random media 

Another example of critical behavior in the presence of weak quenched disorder 
is found for the scaling properties of long flexible polymer chains in random me- 
dia. The question of how linear polymers behave in disordered media is not only 
interesting from a theoretical point of view, but is also relevant for understand- 
ing transport properties of polymer chains in porous media, such as enhanced oil 
recovery, gel electrophoresis, gel permeation chromatography etc. 

It is established that the universal scaling properties of long flexible polymer 
chains in a good solvent are perfectly described within a model of self-avoiding walks 
(SAWs) on a regular lattice Ea. Their scaling properties in the limit of an infinite 
number of steps are governed by the same scaling laws as the joire (undiluted) 
Stanley model at the critical point in the formal to — > limit Ej. In particular, 
for the average square end-to-end distance (-R^) and the number of configurations 
Zn of a SAW with N steps on a regular lattice one finds in the asymptotic limit 
N ^oo: 

(R^) - iV^^ Zn - z^7VT-i (16) 

where and 7 are the universal correlation length and susceptibility exponents 
for the TO = model that only depend on the space dimensionality d and z is a 
non-universal fugacity. For d = 3 the exponents are given in Table 1. 

The problem of SAWs on randomly diluted lattices, which may serve as a model 
of polymer chains i n a porous medium, has been the subject of intensive discussion 
The main question of interest here is: does a small amount of 
quenched structural point-like defects in the medium induce changes to the universal 
properties of the polymer macromolccules? One should be careful in applying the 
Harris criterion "naively" to the SAW problem. Althoueh the critical exponent a of 
a SAW on the d = 3 dimensional pure lattice is positively {a{d = 3) = 0.235±0.003), 
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a weak quenched short-range-correlated disorder does nnt alter the SAW critical 
exponents. This statement haa^been proven by Harris Ell and confirmed later by 
renormalization group results e3. Numerical results for these systems available from 
Monte-Carlo simulations also lead tQ_th£ conclusion that such disorder does not 



affects the scaling behavior of SAWs 

On the other hand, un iversa l properties of SAWs on a percolation cluster obey 
new scaling laws (see, e.g. e3Cj and references therein). 

Another interesting problem is the scaling behavior of polymers in disordered 
media when the defects are correlated or belong to some sponge-like structure. 
In Refs. BEl it was found that the asymptotic behavior of SAWs in long-range- 
correlated disorder of type (|l^) is governed by a set of critical exponents which are 
different from that in the pure case. This recent result waits for its experimental 
verification. 



3.7. Random anisotropy magnets 

Among disordered materials containing magnetic elements there are some inter- 
metallic amorphous alloys, which are characterized not by a random-site disorder, 
but by a disorder in the form of a random anisotropy. Random anisotropy is the most 
important characteristic of these metals and that results from their random close- 
packing structure. Typical examples of random anisotropic materials are amorphous 
rare-earth-transition metal alloys Il3. Their features are taken into account by the 
random anisotropy model (RAM) in ([T^). The model may be extended also to 
describe other intermetallic amoiphous alloys containing rare-earth elements with 
asymmetric charge distribution Eil. 

The original example for the RAM was TbFe2 liS. Some of its structural features 
were used for testing the randomly packed hard spheres model corresponding to 
rare-earth-transition metal alloys. The Tb — Tb correlation function, obtained on 
the basis of generated clusters with spheres of two sizes in Ref. E3 is in qualitative 
agreement with neutron scaitering data on amorphous TbFe2. Using the simple 
charge model, both classical Lj and quantum mechanical calculations O carried out 
on these clusters show the evidence of the anisotropy term in (|l^) as well as that the 
distribution of directions of the local easy axes is effectively isotropic, and indicate 
the apparent absence of correlation between the directions of axes at neighboring 
sites. 

Magnetic properties of the surroundings of rare-earth atoms are neglected within 
this model and thus the RAM does not describe correctly alloys of rare-earth - Fe 
type. Realistic models would have to include both magnetic iron-iron and rare- 
earth - iron interactions. Therefore the RAM is most adequate to describe the 
intermetallic alloys of rare-earth with non-magnetic compounds. 

A considerable number of various experimental data on alloys described by the 
RAM are listed in an exhaustive review by Cochrane et al. lid. As follows from their 
analysis, some of these materials order magnetically at sufficiently low temperature. 
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but the nature of this ordering is not clear. Among possible low-temperature phases 
they discuss ferromagnetic ordering, spin-glass phase, quasi long-range ordering. 
It is clear, that the type of ordering is connected with the alloy intrinsic structure, 
that corresponds within the model to the distribution of anisotropy axes. 

Another interesting application of the RAM is the description of the isotropic- 
nematic phase transition in liquid crystals in mrous media. Monte Carlo simulations 
give the quasi long-range order in this case IZil, while experimental results suggest 
that the bulk nematic phase is replaced by a "glassy" state LS. 



4. Effective Hamiltonians and the renormalization group 

In the theory of critical phenomena it is standard to rely on the methods of renor- 
malization group (RG). The basic conjecture of this approach states that a system 
undergoing a second-order phase transition is scale-invariant at the critical point 
Tc- Near to the critical point the behavior of thermodynamic and correlation func- 
tions are determined only by the correlation length ^ (c.f. Table I) which describes 
the scale of the correlations of the order parameter fluctuations. 

Among the various mathematical schemes of renormalization group the most 
accurate numerical results for critical phenomena are obtained in terms of the field- 
theoretical formulation □. In this formulation the universal critical properties of a 
statistical system are obtained from the large-scale behavior of the corresponding ef- 
fective field-theoretical Hamiltonian. The latter is treated with the well-elaborated 
methods of quantum field theory. The relevant global properties of a microscopic 
model is represented by the structure of the effective Hamiltonian. In particular, 
the effective Hamiltonian that describes the critical behavior of the Stanley model 
(||) reads 

n{$) = J d'':,[l[\S7$f + + (17) 

and is usually referred to as the Hamiltonian of the Landau-Ginzburg- Wilson the- 
ory. Here, is a bare mass proportional to the temperature difference r to the 
critical point, (/) = (/'(R) is a m-component vector field, uq is a bare coupling. Note, 
that the effective Hamiltonian preserves the 0(m) symmetry of the Stanley 
model d). 



4.1. Weakly diluted Stanley model with point-like uncorrelated 
defects 

A weak quenched disorder term can be introduced directly into the effective Hamil- 
tonian (|l^) . As it was mentioned in Section ^ the presence of non-magnetic impu- 
rities in a microscopic model (||) manifests itself in fluctuations of the local tem- 
perature of the phase transition. Introducing ip — ipCR) as the field of local £ritical 
temperature fluctuations, one obtains the effective disordered Hamiltonian 13: 

.2 



d'^a 



{^[|V0? + (A^^ + ^)l0?]+|l^T} 



(18) 
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The Hamiltonian ( p^ ) depends on a number of macroscopic parameters that de- 
scribe the specific configuration of the field tpCR). On the other hand, the observ- 



the random field ip and are 
. In particular, the singular 



ables should not depend on the specific realization o: 
to be averaged over the possible configurations of ■0 
contribution to the free energy of the diluted quenched Stanley model (||) can be 
written in the form of a functional integral: 



Foe D[ij{x)]P[tl;{x)]lnZ[tl;{x)], (19) 



where the configurationally-dependent partition function is the normalizing 

factor of the Gibbs distribution with effective Hamiltonian (|l^); defines a 

probability distribution of the field iIj{x). 

In order to avoid the averaging of the logarithm in the standard step is to 
use the replica trick Lj. This amounts to writing the logarithm in the form of the 
following limit: 

\nZ= lim . (20) 

n^Q n 

While the powers of Z can only be evaluated for integer values of n, analytical 
continuation in n is assumed to perform the limit n — > 0. Introducing n replicas of 



the model (18) and taking that ip obeys a Gaussian distribution: 



P(V') = ^^exp(-V'V4t^') (21) 



with uj'^ being the dispersion parameter, one ends upEa with the effective Hamil- 
tonian: 

2 



— 1 ' — 1 



a— 1 a— 1 




In the limit n ^ the field theory ( |22| ) describes the critical properties of the weakly 
diluted Stanley model. Here, the bare coupling uq is positive (being proportional 
to Uq) whereas the bare coupling vq is proportional to minus the variance of the 
random variable ip and thus it is negative. The last term in ( p2[ ) is present only 
for non-zero dilution: it is directly responsible for the effective interaction between 
replicas due to the presence of impurities 119. 



The model (^2|) is interesting also in the polymer limit m — s- 0. In this case it 
can be interpreted as a model for SAWs in disordered media. Note, that this limit 
is not trivial. As noticed by Kim once the limit m,n ^ has been taken, both 
Uq and vq terms are of the same symmetry, and an effective Hamiltonian with one 
coupling Uq = uq + vq of 0{mn = 0) symmetry results. This leads to the conclusion 
that weak quenched uncorrelated disorder does not change the universal critical 
properties of SAWs. 
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4.2. Weakly diluted Stanley model with long range correlated 
defects 

In the above subsection we obtained an effective Hamiltonian for systems with un- 
correlated point-hke disorder. Now let us consider the case, when the impurities are 
long-range-correlated li^. In particular, let us consider the case, when the inhomo- 
geneities in the system cause fluctuations in the local transition temperature Tc(af), 
characterized by the correlation function g{\x — y\) = {Tc{x)Tc{y)) — {Tc{x))'^, that 
falls off with distance according to the power law (p^. The Fourier-transform g(k) 
of g{x) reads for small k: 

g{k)^vo + wok^-''. (23) 

Note, that in the case of random uncorrelated point-like defects the site-occupation 
correlation function reads: g{\x — y\) ~ 5{x — y), so its Fourier transform obeys: 

5(fc)~«o. (24) 

As above, we apply the replica trick to average the free energy over different 
configurations of quenched disorder and construct the effective Hamiltonian of the 
Stanley model with long-range-correlated disorder OJ: 

n[$) = / d'^x[\Y^ [|v<A"p + + 9 E (ir P)' } + 

E U'xd''yg{\x-v\Wi.^)?\$Hy)?- (25) 

a,/3=l 

Here, the notations are as in formula (p^), and the replica interaction vertex g{x) 



is the correlation function with Fourier image (23). 

Passing to the Fourier image in ( p5| ) and taking into account eq.(p3|), an effective 
Hamiltonian results that contains three bare couplings uo,wo,'Wo. Note, that as far 
as we are interested in the long-wave limit fc — > 0, for a > d the wo-tcrm in ( [2^ ) 
becomes irrelevant and one obtains an effective Hamiltonian of a quenched diluted 
(short-range correlated) Stanley model E2I with two couplings uo,vq. For a < d we 
have, in addition to the momentum-independent couplings, the momentum depen- 
dent one wok°'~'^. Note, that magnetic systems with parallel extended impurities 
are not mentioned in our following description; we refer the interested reader to the 
Ref. for review and numerical results. 

The anahisis of SAWs in media with long-range-correlated disorder uses Kim's 
observation t3 that in the limit m, n — > one may pass to an effective Hamiltonian 
with only two couplings Uq = uq + vq and wq (in what follows below we will keep the 
notation uq for thi s new coupling Uq). In discrete momentum space this effective 



Hamiltonian reads EjEJ: 

n _j 71 4 

k a a kik2k3ki i=l 
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n 2 4 

E \kr'5{Y,h + k)S{Y,h~kX$tM^i$i),ni,n^0. (26) 

a/3 kk-ik2k3ki i—1 i—3 



Here, the 6{k) represent products of Kronecker symbols and the notation 
implies a scalar product. Note that the wo-term contains the interactions between 
replicas and an additional power of an internal momentum. Again, it may be shown 
that for a = d in the limit m,n —> both uq and wq terms are of the same symmetry 
^d one is left with an 0{mn — 0)-vector model with only one coupling (uq + wq) 

4.3. Stanley model with random anisotropy 

In the cases considered so far, point-like defects introduced random scalar variables 
to the effective Hamiltonian ([l7|). Random anisotropy disorder on the other hand 
( ^tI) is different in this respect. The random variable of the disordered effective 
Hamiltonian for the RAM has the form of a vector with random orientation and 
given magnitude. Assuming this magnitude as constant one obtains the replicated 
Hamiltonian as follows 

d'M^EN^"i'+A^orP]+9Ei^^"i'-E^(^^^") • (27) 

I Q = l ■ Q = l Q = l J 

Here, the notations are as in (p^), D is the magnitude of the anisotropy field, z is 
a random unit vector pointing in the direction of the local anisotropy field. 

As it was noted above, to perform the average over the random variables z 
we will use the distributions (|l4| ) and ( |l5|) given in Section ^. For the isotropic 
distribution of the directions of random anisotropy averagingthe axis of all sites 
(|4|) results in an effective Hamiltonian of the following form 



^ f n n 

H(0)= /d'^Ji^[iv0"p+^^irp]+^E 




ia\4 



+ E E«0f^f , (28) 

• a,/3=lij = l J 

with the bare couplings uq > 0, uq > 0, wq < 0. Furthermore, the values of vq and 
Wq are related to D and appropriate cumulants of the distribution function ( p^ ) in 
such a way that their ratio equals wq/vq =^ — m. Note that the symmetry of the 
Uq and Vq terms corresponds to the random site Stanley model (^). However the 
UQ-term has the opposite sign. 

Now, we consider the case of the distribution of z along cubic axes (^5|). After 
performing the averaging over configurations the corresponding effective Hamilto- 
nian reads 111: 
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Q = l / ■ i=l a, 13=1 ■ i=l a=l I 

where the bare coupUngs are uq > Q,vq> 0, wq < 0. The symmetry of the wq terms 
in ( p9[) differs from that in (p8|). Furthermore, the values of wq and vq differ for the 
Hamiltonians (^) and (^9|) due to the cubic distribution ( [l5| ) but their ratio equals 
— m again. The last term in ( p9| ) combines the symmetries of terms with coefficients 
uq and Wq. It does not result from the functional representation of the free energy 
but is generated by further application of the RG transformation. Therefore yo can 
be of either sign. 



4.4. Field-theoretical renormalization group approach 

In this subsection, we give a brief account of the main relations of the field- 
theoretical renormalization group (RG) formalism that can be evaluated in dif- 
ferent variants. We choose the massive field theory scheme with renormalization of 
the one-particle irreducible vertex functions_rg^'^''(fci, .., kL]Pi, --tPn] IJ^t] 
non-zero mass and zero external momenta 113. The one-particle irreducible vertex 
function can be defined as: 

(02(ri) . . . . . . 0(i?jv)>rpf d'^i?i . . . A^RnA^ti ■ ■ ■ dV . (30) 

Here, {Aq} stands for the set of bare couplings of the effective Hamiltonian, {p}, {k} 
are the sets of external momenta, Aq is the cutoff, and the averaging is performed 
with the corresponding effective Hamiltonian, Tie//- 

An intrinsic feature of the evaluation scheme (^) are divergences in the limit 
Aq — !■ oo. In order to remove them and to map the divergent mathematical objects 
to convergent physical quantities one performs a controlled rearrangement of the 
series for the vertex functions. Introducing renormalizing factors for the fields Z^, 
Z^2, the renormalized vertex functions F^'''^'' are expressed in terms of the bare 
vertex functions as follows: 

rif '^^{fc}; M; m'; {A}) = z^.z^/'ri'^'''\{k}; M; {Ao}), (3i) 

where {A} are the renormalized mass and couplings. 

The regularization scheme suggested in Eq. (^Tj) is not unique. To define the 
regularization scheme completely one imposes renormalization conditions for the 
renormalized vertex functions. In the massive scheme the renormalization conditions 
read 10: 

Tf^\k,-k■,p^{x})\^^o^^^^ 

-±rf'\k,~k;p^{X})U^o = h 
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T^^^^\p;k,-k;^,^,{X})\k=p=o = 1- (32) 

In the case of long-range-correlated disorder we have another global parameter a 
along with d, and for the renormalization of the coupling w one imposes 0: 

Tf;t\{k};i^Mmk=o^f^'-''w (33) 

with the renormalized coupling w. The scaling properties of the system asymptoti- 



cally close 



o the critical point are expressed by the homogeneous Callan-Symanzik 



equations for r'-^'^\{p}; {k}; n'^, {A}): 

(|-L)7^({A}KL70^({A})}^^^■^)(M;{fc};^^{A}) = O, (34) 
with the coefficients: 

a ^^\^\ i ^^"P i - dZ^2 

- 9h^'{^°>-^'" 9h^'{^°>^''°' " ~ah^'{^°>-^° ^^^> 

where = Z^'iZ^f,. 

The fixed points {A*} of the RG transformation are given by the solution of the 
system of equation: 

/3a,({A*}) = 0, z = 1,2,.... (36) 

The stable fixed point, corresponding to the critical point of the system, is defined 
as the fixed point where the stability matrix: 

possesses eigenvalues {wi} with positive real parts. In this point, the function 7^ 
determines the value of the pair correlation function critical exponent 77: 

'/ = 70({A*}), (38) 

and the correlation length critical exponent u is determined as: 

7.-i=2-74{A*})-70^({A*}). (39) 

All other critical exponents may be obtained from familiar scaling laws. For exam- 
ple, the susceptibility exponent 7 is given by: 

7 = 1^(2-7?). (40) 

The /9- and 7-functions are calculated perturbatively as series in the (several) 
couplings Xi. The order of the expansion corresponds to the number of loops in 
the diagrammatic Feynman representation of the vertex functions (|3^). However, 
due to the (asymptotic) divergence of the RG functions series, one cannot directly 
derive the reliable physical information from these expressions. The methods that 
permit to cope with the task will be described in the next section. 
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5. Resummation of asymptotic series 

The expansions of RG functions in powers of couplings (weak coupling expansions) 
are known for a number of models with high accuracy. However, the very idea of 
perturbation theory is that the result is successively approximated by accounting 
higher order contributions. But the method seems doomed to failure if the series are 
divergent. This is the case in the field-theoretical renormalization group approach. 
Moreover, here the series are characterized by a factorial growth of the coefficients 
implying a zero radius of convergence. A simple method to manage this problem 
is optimally truncating the series, i.e. accounting only its first several terms that 
do not show the divergence. A self-consistent way to take into account the higher 



resummation: Borel 
Resumming a series 



order contributions requires the application of special ways 
resummation accompanied by certain additional procedures r ^ 
one should note that the accuracy of the resummed results depends on the proce- 
dure of resummation. Thus the analysis and appropriate selection of an adapted 
resummation procedure is crucial to obtain reliable results. 

However, there remains the principal question about the Borel summability of 
the perturbation theory series for a given model. Up to now this has been proved 
only for the (j)'^ theory with one coupling Ej. For years, the field theoretical RG 
series for models with several couplings were analyzed as if they are asymptotically 
divergent without a proof of this property. Moreover, there exists strong evic 
possible Borel non-summability of the series obtained for disordered models li 
In this section, we will give the definition of an asymptotic series and provide the 
methods which currently are used to extract convergent data. 



1/(5) -Jl^fef'"'! - 



5.1. An asymptotic series 

For a function f{g) given in the form of a power series 

OG 

fig) = Y,"-kg\ 

k=0 

the basic definition of its asymptotic nature is: 

N 

-tfey I :i "fe+i; 
fc=o 

for fixed N and small g. 

In general, the nature of an asymptotic series is such that the correct sum is 
uniformly approached until, after an optimum number A^o of terms the subsequent 
terms drive thejsartial sums away from the correct result for the total sum, leading 
to divergence H. 

For models with a simple 0(m)-symmetry the factorial growth of the RG. func- 
tions coefficients for large orders of perturbation theory is well established ES: 

ak^ k\{~afk''c[l + Oiy)], k ^ 00 (41) 
k 
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the quantities a and b have been calculated in Ref. Ezl, the values of constant c 
in Ref. @. The estimate (^ ) is a fundamental one since it shows that the initial 
functions can be reconstituted from their asymptotic expansions by the application 
of the integral Borel transformation in different modifications. As it was already 
noted above the asymptotic nature of the divergence of the RG functions for the 
disordered models still has not been proved. 

Below we will describe the methods of resummation which in the field-theoreti- 
cal RG scheme are applied to study the critical behavior of the Stanley model ( p7[ ) 



as weU as the disordered models (M), (M), (E§ 



5.2. Resummation of asymptotic series with one coupling 

The idea of resummation which can be applied to an asymptotic series consists 
in changing the order of summation0. In the case of one variable it means that 
starting from a power series 

L 

S{x) = lim y^aix\ (42) 

i=l 

one performs an identical transformation 

S (x) EE hm^^oo ir dteM^m ^ 

hmi^oolliniA^oo /o" rfiEti ^ exp(-t)]}, (43) 
and redefines the sum by the following expression: 

S' {x) = limA^oo{liniL^oo ^^Ef^i ^^^^ exp(-t)]} = 

j^dteM-t)i:t,^-^, (44) 

where X^i^i CLi{xty /i\ is called the Borel - image of o-i^^^ ■ This trick is natural 
in the sense that for the case of a convergent series one finds S' — S within the 
radius of convergence. 

The above mentioned procedure in the case of a single variable is known as 
the Borel resummation technique and in different modifications it is widely used 
for evaluating asymptotic series. Unfortunately this technique cannot be applied 
to our case immediately because only the partial sum of the truncated series is 
known. To overcome this obstacle one represents the Borel-image of the initial sum 
in the form of a rational approximant and in this way reconstitutes the general 
form of the series. This technique that involves a rational approximation_and the 
Borel transformation is known as the Pade-Borel resummation technique l2I and is 
performed as follows: 

• the Borel-image of the initial sum is constructed: 

E«^-'^E^; (45) 
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• the Borel-image is extrapolated by a rational approximant 

[M/N] = [M/N] {xt); (46) 

here, [M/N] stands for the quotient of two polynomials in xt; M is the 
order of the numerator and N is that of the denominator; 

• the resummed function is then obtained by the integral: 

S'''"'{x)= dtexp{~t)[M/N]{xt). (47) 



Jo 

The Pade-Borel procedure can be generalized by introducing an additional fit pa- 
rameter p to the Borel transformation. Substituting the factorial i\ by the Euler 
gamma- function T{i + p + 1) and inserting an additional factor into the integral 
(|47|), one defines the Pade-Borel-Leroy resummation procedure. 

The Pade-Borel as well as Pade-Borel-Leroy resummation procedures are ap- 
plicable to RG functions of one coupling constant such as those of model (pT|). 
For models with more complicated symmetries and several couplings one uses the 
generalizations of the Pade-Borel method as explained in the following subsection. 

5.3. The case of several couplings 



In order to suit the resummation procedure (45)-(^7|) for functions that depend 
on several variables one should change the first_step (|45|): for example, for the 
two- variable case one defines the Borel image by Eil: 



The rational approximation ( ^ ) can be performed then either in the dummy vari- 
able t or in the variables x, y. This defines two approaches that one can use to resum 
the functions of several variables. In the first case the resurnmation procedure is 
referred to as the Pade-BoreLxfisummation for resolvent series e3. The application of 
the Chisholm approximants Ed which are the generalization of Pade-approximants 
to the many-variable case is necessary in the second case. A Chisholm approximant 
can be defined as a ratio of two polynomials both in variables x and ?/, of degree M 
and N such that the first terms of its expansion are equal to those of the function 
which is approximated. Again, the resummation is performed in eq. ( p7[ ) replacing 
the Pade with the Chisholm approximant. This method will be referred below as 
the Chisholm-Borel resummation. 

It is obvious that an initial sum can be resummed in different ways. Apart 
from the Leroy fit parameter p mentioned above some arbitrariness arises from the 
different types of rational approximants one may construct. For instance, within 
the two-loop approximation the method of Pade-Borel resummation of a resolvent 
series can be done using either the [0/2] or the [1/1] approximants. The Chisholm- 
Borel approximation implies even more arbitrariness and demands a very precise 
analysis of the approximants to be chosen. 
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We want to stress here again that in contrast to the single variable case the 
validity of neither the Pade-Borel nor the Chisholni-Borel resummations have been 
proven rigorously for any cases with two or more variables. The success of their 
application can solely be judged by the consistency of the results with both exper- 
imental and Monte-Carlo data. 

5.4. The method of subsequent resummation 

In the resummation schemes described above for series in several variables the 
variables are treated as "equal in rights" . This is not the case in the recently pro- 
posed method of so-called subsequent resummation. developed in the context of the 
d = 0-dimensional disordered Ising model in Ref. £3. 

Starting from the RG function f{u, v) in the form of a series in powers of u, v 
(for the model (|2^), u is the original coupling of the undiluted system and v is the 
variance of the quenched disorder) 

/Kz;) = ^^Cfc,„uV, (49) 

n k 

one rewrites it as a series in the variable v: 

/(«,,;) = ^A„ («)«", (50) 

n 

with coefhcients that are in turn series in the variable u and are to be resummed 
in advance by any one- variable method: 

A„(u) = ^Cfc,„u^ (51) 

k 

The main result of Ref. S is that the expansions of the coefhcients ( ^ ) and the 
resulting series ( |50| ) at fixed u are Borel summable (for a d = dimensional system, 
however). This gives a hint to analyze the RG expansions of the d = 3 disordered 
models by first performing a Pade-Borel resummations of the corresponding se- 
ries for the coefficients in one coupling and then, using the computed coefficients, 
resumming the series in the other subsequent coupling(s). 

6. Results of the renormalization group analysis 

Now, having at hand the powerful methods of the field-theoretical renormalization 
group refined by resummation of the (asymptotic) series we can check to what 
results does it lead when applied to systems with weak quenched disorder. Again, 
we start from the uncorrelated point-like disorder, then pass to the systems with 
long-range correlated disorder, and finally we analyze random-anisotropy systems. 
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6.1. 



Weakly diluted Stanley model with point-like uncorrelated 
defects 



Referring to the universal critical properties of the model (22) at d = 3, as implied 
by the Harris criterion one observes new critical behavior only for the case m ~ 1, 
i.e. for the weakly diluted quenched Ising model. 



Following the procedure sketched in the subsection 4.4 the expressions for the 
RG functions of the model are obtained as series in the renormalized couplings 
u and V. Each order of the perturbation theory corresponds to a given number 
of integrations in momentum space and, equivalently, to the number of loops in 
Feynman diagrams in the diagrammatic representation for the vertex functions 

The massive RG functions of the weakly diluted quenched Ising model are known 
to high orders of the perturbation theory. They result from more than twenty years 
of laborious studies. Within the three-loop accuracy the functions were first re- 
ported in Ref. Ej. They yet contained errors which were corrected in Ref. l3, er- 
roneously again. The final three-loop expressions were only published in Ref. l3. 
Subsequently, four-loop series were obtained in Ref. tH. Recently five-loop l3 and 
record six-loop @ expansions became available. Written in the two-loop approxi- 
mation the functions read E£j: 



(3u{u,v) 
702 {u,v) 



1 



-(4- 

-[I2(n-- 



3 

— V - 

2 



l2\UV 



8 1 

4L V 1 



(52) 



1 — V u - 

3 



1 



-)+i2\vu+— [3(n - -)+i2\'' 



-2(4 
(4-^){i 



16 



«2 



-V - 12 
4 



27 



1 

— 1 
12 



32 



(53) 
(54) 

(n-i)}. (55) 



Here, we use the conventional normalization for the couplings u, v such that the 
coefficients of u and v in the one-loop contribution to /?„ and [3^ equal unity. 

The expressions (|5^)-(|55|) besides the couplings w, v depend on space dimension 
d. In particular, d enters the two-loop integrals ii(d), i2{d). The two-loop contribu- 
tions to the RG functions of the three-dimensional weakly diluted quenched Ising 
model can be cjM^ined from the expressions (p2|)-(55) by substituting ii(3) = 2/3, 
i2(3) = -2/27 



The explicit depe ndence of the functions (|52|)-(|55D on the space dimension d was 



exploited in Ref. Elj. Expanding the values of the loop- integrals ii, 12 in e = 4 — d, 
the functions were obtained in the form of an e-expansion E3. Classically, this result 
follows the familiar e-expansion IihI of the critical exponents. In order to obtain the 
expansions, one solves the fixed point equations (|3^) with e as a small parameter 
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and then re-expands in s the critical exponents after substitution of the coordinates 
u*{e), v*{e) into the 7- functions. 

This form of analysis of the six-loop massive RG functions encounters hardships 
since the e-expansions of the loop integrals are not known in high four-, five- and 
six-loop approximations. This obstacle can be overcome by anpb-ing an alternative 



dimensional regularization scheme with minimal subtraction tll3. In this approach 



the RG functions were obtained in the three-loop cl!3 approximation and can be 



recons titu ted within the five-loop accuracy from the RG functions of the cubic 
model 1123. The degeneracy of the one-loop R G func tions lead to the result, that 

-^,106 



instead of expanding in e one is to expand in ^/e 



This results in y^i^pansions 



for the critical exponents E!!Z1 and the stability matrix {p7h eigenvalues 



108 



V = 0.5^-0.08411582el/2_o,ol663203e-h0.04775351e3/2_^0.27258431e^ 
7] = -0.00943396e-h0.03494350e^/2-0.04486498e2-H0.0215732l£5/2 , 
7 = l+0.16823164el/2-0.02854708e-h0.07882881e3/^-h0.56450490e^ 
Loi ^ 2 £+3.70401119 e^/^ll. 30873837 e^ 



UJ2 = 0.67292659 £^/2-1.92550909 £-0.57252518 £ 



3/2_ 



13.93125952 £^ 



(56) 
(57) 
(58) 
(59) 
(60) 



In the above expressions the two-loop results for the exponents were obtained in 



and 



LIO 



Ref. L3, the three-loop results were presented independently in Refs. 

The crucial property of the •y/£-expansions (|56|)-(|60|) is their formal character. 
Though they allow t o pr edict qualitatively new critical behavior in the weakly 
diluted Ising model L3£23 they seem to be of no use for any quantitative analysis. 
Naively adding the successive perturbational contributions in the y^-expansion for 
the model stability matrix eigenvalues (|59|), ( |60| ) one observes that already in the 
three-loop approximation (~ e) wojiecomes negative and therefore no stable fixed 
point exists in stri ct \/£-expansion £23. Even resummation procedures do not change 
this picture E3lli] . 

The way to cope with the bad convergence properties of the -y/E-expansion can 
be found in the fixed-dimension approach which does not imply expanding the loop 
integrals in e. This method yields most numerical results for the universal charac- 
teristics of the diluted Ising model at criticality. The numerical solution of the fixed 
point equations (^) the resummed /3-functions ( |5^ ) and the resummation of the 7- 
functions ( |55| ) in the stable fixed point shows the validity of the method. The study 
of the massive /3-functions of the diluted Ising model resummed in this way revealed 

that starting from the two-loop a.pproximatioyi |.he ra^ndom^|ixed P p'n t is sta ble and 

i,U^, 115|ll6 ,ll^ 



is present in all orders of perturbation theory 



on the three-dimensional values of the loop-integrals llHil, the RG functions of the 



Based 



model are calculated and resummed 1113 at present in the six-loop approximation. 



Let us analyze the fixed point structure of the /3-functions for the random Ising 
model at d = 3. Within the one-loop approximation (i.e. including only linear terms 
in u, V in the curly brackets in (p^, (^) at d = 3 one finds three fixed points 
(see Fig. ||): the Gaussian fixed point G u* = v* = 0, the pure Ising fixed point I 
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Fig. 6. Qualitative RG flow for the weakly diluted quenched Ising model. The Gaussian fixed 
point G is stable for d> 4, stable fixed point P can not be reached from the initial coupling values 
M > 0, i; < 0. The unphysical regions for the random Ising model are shown in grey tone. The 
pure Ising fixed point I is unstable whereas random fixed point R is both stable and accessible 
(stable fixed points are shown by boxes). 



u* ^ 0, V* — and the "polymer" fixed point P u* — 0, v* 0. It is straightforward 
to check that the fixed points G and I are unstable whereas the fixed point P is 
stable but as far as > it is inaccessible for the initial values of couplings of 
our model. Within the one-loop approximation one does not encounter the fixed 
point R with both non-zero coordinates u* ^ Q, v* ^ 0: this happens becauaaJJi 
system of equations for the fixed points is degenerate on the one-loop level L§E2^ 
As we will see, this fixed point appears in the next, two-loop approximation if the 
resummation is applied. 

Applying the Chisholm-Borel resummation procedure of the subsection 5.3 one 
encounters the random fixed point R of the model in the two-loop approximation 
(see Fig. The stability analysis shows that this fixed point is stable proving the 
crossover to a new critical regime under dilution. In Figs. ^ ^ we show the curves 
I3u{u,v) = 0, f3y(u^v) = in the u — v plane. The intersections of these curves 
(i.e. simultaneous zeros of both /3-functions) correspond to the fixed points. The 
"naive" analysis of the /3-functions, without applying any resummation procedure 
leads to the curves, which are shown on the left-hand side of the Figs. |^, ||. Without 
resummation only in the three-loop approximation one finds a stable random fixed 
point u* ^ 0,v* ^ 0. However the fixed point then disappears in the four-loop 
approximation. A completely different picture is observed when the resummation 
procedure is applied (right hand columns of the figures). In the region of interest 
for the values of the couplings the topology of the fixed point picture remains stable 
for increasing orders of the approximation from the two-loop to the four-loop level. 
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Fig. 7. The lines of zeros of the diluted Ising model massive /3-functions in the one- and two-loop 
approximations of the perturbation theory, left-hand column: non-resummed, right-hand column: 
resummed by the Chisholm-Borel method. Circles correspond to I3u = 0, thick lines depict /3„ = 0. 
Thin solid and dashed lines show the zeros of the analytically continued functions /3u and 
respectively. One can see the appearance of the mixed fixed point m > 0, ''r^O in the two— loop 
approximation for the resummed /3-functions. The figure is taken from Ref. 



The Figs. 0, ^ show also the historical evolution of the study of the weakly diluted 
quenched Ising model. Already in the two-loop level the RG functions resummed 
by the Chisholm-Borel procedure 1113 revealed that no difficulties connected to the 
degeneracy of the /3-functions are encountered. Critical exponents extracted from 
the resummed 7-functions values in the fixed point were found to be clearly larger 



above, in the three-loop level 
yields fixed point coordinates 



than those of the pure model (see Table ^). As 
the straightforward analysis of the /3-functions 
and critical exponents without resummation, but the accuracy obtained did not 
allow to estimate, for instance the heat capacity critical behavior. Numerical values 
of the random Ising model critical exponents as obtained by the RG approach are 
given in T able ^ with reference to the renornialization schemes and resummation 
procedures 1113. We note here a breakthrough that occurred during several months 
of 1999-2000 wh£n the perturbation [hjnry expansions were extended from the 4th 
lZI through 5th lB to the sixth order 



Already the 5th loop order Padc-Borcl resummation meets some difficulties. 
Moreover, the extension to the six-loop order revealed a wide gap between the 
five- and six-loop fixed point coordinates tlJ and subsequently an inconsistency 
of thesix-loop values of critical exponents compared with the five-loop results of 
Refs O. One also encounters a similar behavior of the resummed series when the 
d = 3 minimal subtraction RG technique is applied to the random Ising model: 
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Fig. 8. The lines of zeros of the diluted Ising model massive /3-functions in the three- and four- 
loop approximations of the perturbation theory, left-hand column: non-resummed, right-hand 
column: resummed by the Chisholm-Borel method. The notations are the same as in the figure 
^ Close to the mixed fixed point the behaviour of the resummed functions remains alike with 
the increasing order of a]pp*pximation. This is not the case for the non-resummed functions. The 
figure is taken from Rof. Eij. 



the random fixed point is present in the resummed two- three- and_EQiir- loop 



111 



This 



predicted already in a study of the randomly 
^. The perturbation theory series of this toy- 



approximations whereas it disappears in the five-loop approximation L3t 
lead to the conjecture about the possible Borel-non-summability of the series. 

The fact that the weak coupling expansion might not be of asymptotic nature 
for the random Ising model was 
diluted model in zero dimensions 
model appeared to be Borel non-summable. Moreover, this property was shown to 
be a direct consequence of the appearance of Griffiths-like singularities l5 that are 
caused by the zeroes of the partition function of the pure system. 

However, as noted above, in another recent study of the d — dimensional 
random Ising model it was shown analytically that the perturbative expansions 
for the free energy are Borel summable provided the summation is carried out 
subsequently, as described aboveiii subsection 5.4. The application of this procedure 
to the case d = 3 reconstituted tia earlier data for the asymptotic critical exponents 
of the model in massive scheme and allowed to restore the presence of the xaiidom 
fixed point in the d = 3 minimal subtraction scheme on the five loop level IIHI . 

Alternatively, the d = 3 random Ising model was treated non-perturbatively 



first by means of the Golner-Riedel scaling field approach LLHil, a.nd very recently by 



an RG approach based on the concept of effective average action 



122 



February 1, 2008 18:40 WSPC/Guidelines lecture 



32 Yu. Holovatch, V. Blavats'ka, M. Dudka, C. von Ferber, R. Folk, and T. Yavors'kii 



Table 3. The theoretical values for the asymptotic critical exponents of the weakly diluted Ising 
model. rtLA denotes the nth order in loopwise approximation within the massive ('mass') and 3d 
minimal subtraction ('MS') schemes of the field-theoretical renormalization group approach. The 
resummation procedures are given in the following notations: ChBr - Chisholm-Borel; PdBr - 
Pade-Borel; AW - e algorithm of Wynn, CM — Borel transformation with conformal mapping. SF 
stands for Golner-Riedel scaling field method, EEA denotes nonperturbative RG approach based 
on the concept of effective average action, superscript at the correction-to-scaling exponent ui 
indicates that the real part of the corresponding complex number is shown. 



Ref. 


RG 
scheme 


Order 


Resum- 
mation 


u 




7 


UJ 


Sokolov et 
1981, Ref. 






mass 


3LA 


No 




0.009 


1.31 




Newman et 
1982, Ref. 


121 




SF 


No 




0.70 


0.015 


1.39 


0.41 


Jug, 1| 
Ref. 


112 


i, 


mass 


2LA 


ChBr 


0.678 


0.031 


1.336 


0.450'= 


Mayer et : 
1984, Ref. 


113 




mass 


2LA 
3LA 


ChBr 
ChBr 




0.031 
0.022 


1.337 
1.325 




Mayer et e 
1989, Ref. 


'h 




mass 


4LA 


ChBr 


0.670 


0.034 


1.326 




Mayer, 
Ref. 


114 




mass 


4LA 
4LA 


AW 
PdBr 


0.6680 
0.6714 




1.318 
1.321 




Shpot, 
Ref. 




i9, 


mass 


3LA 


ChBr 


0.671 


0.021 


1.328 


0.359 


Janssen et 
1995, Ref. 


104 




MS, 3d 


3LA 


PdBr 


0.666 




1.313 


0.366 


Holovatch e 
1997, Ref. 


116 




mass 


3LA 


ChBr 


0.671 


0.019 


1.328 


0.376 


Folk et air;— 1 
1998, Ref. Eiil 


MS, 3d 


2LA 
3LA 
4LA 


ChBr 


0.665 
0.654 
0.675 


0.032 
0.022 
0.049 


1.308 
1.293 
1.318 


0.162 
0.430 
0.390^= 


Folk et a 
1999, Ref. 


108 


1 


MS, 3d 
mass 


4LA 
4LA 


ChBr 
ChBr 








0.39(4)'= 
0.372(5) 


Pakhnin et 
2000, Ref. 


)8 




mass 


SLA 


PdBr 


0.671(5) 


0.025(10) 


1.325(3) 


0.32(6) 


Varnashc 
2000 Ref. 


HE 




mass 


4LA 


PdBr 
PdBr 


0.681(12) 
0.672(4) 


0.040(11) 
0.034(10) 


1.336(20) 
1.323(10) 


0.330 


Pelissetto et 
2000, Ref. 


a' 
118 




mass 


6LA 


PdBr- 

CM 
PdBr- 
PdBr 


0.678(10) 
0.668(6) 


0.030(3) 
0.0327(19) 


1.330(17) 
1.313(14) 


0.25(10) 
0.25(10) 


Tissier et 
2001, Ref. 


m 




EEA 


No 




0.67 


0.05 


1.306 





6.2. Systems with long-range-correlated disorder 
6.2.1. Magnets with long-range correlated disorder 

Now, let us discuss the results obtained for the m-vector model with long-range- 
correlated disorder described bv the effective Hamiltonian (25). A one- loop ap- 
proximation was given in Ref. Ej using an expansion ine = 4 — rf, (5 = 4 — a. A 
new long-range-correlated fixed point was found with a correlation-length exponent 
v = 2/a and it was argued, that this new scaling relation is exact and also holds in 
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higher order approximations. This result was questioned recently in Refs. 0, where 
the static and dynamic properties of 3d systems with long-range-correlated disor- 
der were studied by means of the massive field-theoretical RG approach in a 2-loop 
approximation, for different fixed values of the correlation parameter, 2 < a < 3. 
The (3 and 7 functions in the two-loop approximation were calculated as an expan- 
sion series in renormalized vertices u, v and w. The Pade-Borel resummation was 
used for their analysis. The study revealed the existence of the stable random fixed 
point with u* ^ 0, v* ^ 0, w* ^ in the whole region of the parameter a. The 
comparison of the obtained exponents values at various jLand ratio 2/ a shows the 
violation of the supposed exact the relation ~ 2/ a of ref.o. The authors concluded 
that the revealed difference is caused by the use of a more accurate field-theoretical 
description using higher order approximations for the 3d system directly together 
with methods of series summation. Recently, the correction-to-scaling exponent for 
magnets with long-range-correlated defects was obtained using the same approach 



6.2.2. Polymers in media with long-range- correlated quenched disorder 

The effective Hamiltonian (|2^) is the starting point to study the polymer limit 
m ^ of the weakly diluted Stanley model with long-range-correlated disorder. As 
we explained in the subsection 3.6 it may serve as a model for polymers in porous 
media. Again, imposing the renormalization conditions of the massive scheme (^2]), 
( |3^ ) the one-loop approximation leads to the following expressions for the RG 
functions 113: 





4 2/ 








u - - u^h 


— S2uw 





i2S-e)-w^h, 



1^^ = 



2 2/ 
-WI2 
3 



e- [wuh - w'^h] , 



'^S 3" 70 = 0-^ ^4- 



(61) 

(62) 
(63) 



Here, li are one-loop integrals that depend on the space dimension d and the cor- 
relation parameter a: 



h = 



(g2 + 1)2 ■ 



h - 



h 



(g2 



d 

5fc2 



1)"2' = 

[q + kY + l 



d(fg2(a-rf) 

(g2 + 1)2 ■ 



(64) 



Note that in contrast to the usual 0^ theory the 7^ function in Eq. (63) is nonzero 
already at the one-loop level. This is due to the fc-dependence of the integral 1^ in 
Eq. (jel). 

There are essentially two ways to proceed in order to obtain the qualitative 
characteristics of the critical behavior of the model. The first is to consider the 
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polynomials in Eqs. (|6l|), ( |6^ ) for fixed a,d and look for the solution of the fixed 
point equations. It is easy to check that these one- loop equations do not have any 
stable accessible fixed points for d < 4. The second scheme is to evaluate these 
equations in adouble expansion in e = 4 — d and (5 = 4 — a as proposed by Weinrib 
and Halperin EJ. Substituting the loop integrals in Eqs. (|6l|)-(|63|) by their expansion 
in e = 4 — d and 6 ~ A — a , one obtains the 3 fixed points given in the table |[ The 

jible 4. Fixed points and stability matrix eigenvalues in the first order of the e, <5 - expansion 



Fixed Point 


u* 


w* 




Gaussian (G) 








-e -5 


Pure SAW (P) 


e 





e e/2-5 


Long-range (LR) 




S{s-2S) 


^{£-4<5± V£2-4£<5-(-852} 



following conclusions may be drawn from these first order results: Three distinct 
accessible fixed points are found to be stable in different regions of the a, d-plane: 
the Gaussian (G) fixed point, the pure (P) SAW fixed point and the long-range 
(LR) disorder SAW fixed point. The corresponding regions in the a, c?-plane are 
marked by I, II and III in Fig. ||. In the region IV no stable fixed point is accessible. 

For the correlation length critical exponent of the SAW, one finds distinct values 
^'puic for the pure fixed point and i^lr for the long-range fixed point. Taking into 
account that the accessible values of the couplings are w > 0, w > 0, one finds that 
the long-range stable fixed point is accessible only for J < e < 25, or d < a < 2+d/2, 



a region where power counting in Eq. (26) shows that the disorder is irrelevant. 



(65) 



In this sense the region III for the stability of the LR fixed point is unphysical. 
Formally, the first order results for d < 4 read: 

i^p^c - l/2 + e/16, 5 < e/2, 
i^LR = 1/2 + (5/8, e/2<S< e. 
Thus, in this linear approximation the asymptotic behavior of polymers is governed 
by a distinct exponent i^lr in the region III of the parameter plane a, d. However, 
the region where the LR fixed point is found appears to be unphysical. 

Something similar happens if the e, (5-expansion is applied tojitudy models of 
m-vector magnets with long-range-correlated quenched disorder lUS: also in the 
case of magnets, as well as for polymers the first order e, 5-expansion leads to a 
controversial phase diagram. In order to obtain a clear picture and more reliable 
information, one should proceed to higher order calculations. 

To investigate the peculiarities of the critical behavior in the 2-loop approxima- 
tion one may make use of the m ^ limit of the appropriate m-vector mod el ( ^5|) , 
investigated recently 0. Starting from the two-loop expressions of Ref. for the 
RG functions of the Stanley model with long-range-correlated disorder and mak- 
ing use of the symmetry arguments 00 in the polymer limit rn = one gets the 
following expressions for the d = 3 RG functions of the model described by the 
effective Hamiltonian (kG 
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Fig. 9. The critical behavior of a polymer in a medium with long-range-correl£|te|d disorder 
in different regions of the d, a-plane as predicted by the first order e, (5— expansion LJ- Region I 
corresponds to the Gaussian random walk behavior, in the region II scaling behavior is the same 
as in the medium without disorder, in region III the "long-range" fixed point LR is stable and 
the scaling laws for polymers are altered, in region IV no accessible stable fixed points appear; 
this may be interpreted as the collapse of the chain. 



95 1 

(3u{u,w) ^ ~u + - (3/i(a) - f2{a))uw - — + -62(0)7/^^ - 

zib o 

{bM - \h{a))uw'' + h{a)w^ + h{a)w\ (66) 

/3^(w, w) = -{A - a)w - (/i(a) - h[a))w^ + "f" + bw{a)w^ - 
23 1 

— u^w + -6i2(a)uu;^, (67) 
21d 4 

70(m, w) = ^/2(a)w + Y^"^ + ci{a)w^ - ^C2{a)uw, (68) 

%2{u,w) = - - Y^'u^ - C3{a)w'^ + ^C4{a)uw. (69) 

Here, the coefRcients fi{a) are expressed in terms of the one- loop integrals in Eq. 

f ) , hi (a) and Ci (a) originate from the two- loop integrals and are tabulated in Ref. 
for d — 3 and different values of the parameter a in the range 2 < a < 3. The 
series are normalized by a standard change of variables u — > -^Ii, so 
that the coefficients of the terms u, in /3„ become 1 in modulus. 

As in the former cases discussed in subsections ^j^, |6.2.1 the Question about the 
summability of the series in Eqs. (^^ - (^9|) is open. In Ref. Ej various kinds of 
resummation techniques have been applied in order to obtain reliable quantitative 
results for the system under consideration and to check the stability of these results. 
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1 2 3 4 5 6 




1 2 3 4 5 6 

U 



Fig. 10. The lines of zeroes of the 3d /3-functions (pq), ( |67[ ) resummed by the Chisholm-Borel 
method at a = 2.9. The dashed line corresponds to f3u = 0, the solid lines depict f^n, = 0. 
The intersections of the dashed and solid lines give three fixed points shown by filled circles at 
u* = 0,w* = (G), u* = 1.63,10* = (P), and u* = 4.13,10* = 1.47 (LR). The fixed point LR 
is stable. 



First, a simple two- variable Chisholm-Borel resummation technique was employed, 
which turns out to be the most effective one for the given problem. In addition to 
the familiar fixed points G and P describing Gaussian chains and polymers (SAWs) 
on regular lattices, the stable long-range fixed point LR for polymers in long-range- 
correlated disorder was found. In Fig. |l^ the lines of zeroes of the resummed f3- 
functions (|66|), (67) at a = 2.9 in the u, w-plane in the region of interest are depicted. 
The intersections of these curves correspond to the fixed points. The corresponding 
values of the stable fixed point coordinates and the stability matrix eigenvalues for 
different values of the correlation parameter a < 3 are given in Table |^. 

To calculate the critical exponents the same resummation technique was applied. 
The numerical values for 7 and rj are listed in Table || for a = 2.3, . . . , 2.9. Note, 
that for a = 3, which corresponds to short-range-correlated point-like defects, the 
interactions u and w become of the same symmetry, so one can pass to one coupling 
(u — w) and reproduce the well-known values of the critical exponents for the pure 
SAW model. The numerical values corresponding to those listed in Table |^ in this 
case read: u* — 1.63, — 0.59, 7 = 1.17, f] = 0.02, ui = 0.64. As departing from the 
value a = 3 downward to 2 one notices a major increase of the value of the coupling 
u, so the results are more reliable for a close to 3. At some value a = a„iarg the LR 
fixed point becomes unstable. 

To verify the results, the method of subsequent resummation (see subsection 
5.4 ) was applied to the RG functions (66), (67). Here, the summation was carried 
out first in the coupling u and subsequently in w. Again, the presence of stable 
fixed point LR for flmarg < a < d was confirmed. 

The obtained results may be summarized and interpreted as follows: 
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Table 5. Stable fixed point of the 3d 2-loop /3-functions, resummed by Chisholm-Borel method, 
the corresponding critical exponents and the stability matrix eigenvalues at various values of a. 



a 


u* 


UI* 


V 


7 


V 


^1,2 


2.9 


4.13 


1.47 


0.64 


1.25 


0.04 


0.25 ± 0.62 i 


2.8 


4.73 


1.68 


0.64 


1.26 


0.04 


0.22 ± 0.76 i 


2.7 


5.31 


1.81 


0.65 


1.28 


0.03 


0.18 ± 0.89 i 


2.6 


5.89 


1.87 


0.66 


1.29 


0.03 


0.15 ± 0.99 i 


2.5 


6.48 


1.89 


0.66 


1.31 


0.02 


0.11 ± 1.09 i 


2.4 


7.10 


1.87 


0.67 


1.33 


0.01 


0.07 ± 1.18 i 


2.3 


7.76 


1.84 


0.68 


1.36 


0.01 


0.03 ± 1.26 i 



(i) A new stable fixed point (LR) for polymers in long-range-correlated disorder 
is found for d = 3, a < d, leading to critical exponents that are different from those 
of the pure model; 

(ii) There is a marginal value amaig for the parameter a, below which the stable 
fixed point is absent, indicating a chain collapse of the polymer for disorder that is 
stronger correlated. 

(iii) The critical exponent v increases with decreasing parameter a, like in the 
Weinrib and Halperin case. But the relation v = 2la does not hold. Physically 
this means that in weak long-range-correlated disorder (a > flmarg) the polymer 
coil swells with increasing correlation of the disorder. The self avoiding path of the 
polymer has to take larger deviations to avoid the defects of the medium. 



6.3. Critical behaviour of Stanley model with random anisotropy 



The RG results discussed so far in the subsections 3.1, 6.2 concerned random- 



temperature like disorder. Another type of disorder that we consider here is the 
disorder induced by a single-ion anisotropy term in the RAM Hamiltonian (p^. 
According to the chosen distributions of the anisotropy axis (|4|) and ( [l5| ) our 
consideration falls into two parts. First we analyze the case of the RAM with 
isotropic distributions of anisotropy axis and then we discuss what peculiarities in 
the critical behavior of the RAM might appear if the random-anisotropy axis is 
taken to point only along the edges of m-dimensional hypercube (we call this the 
" cubic" distribution) . 



6.3.1. An isotropic distribution of the random anisotropy axis 



Applying the renormalization conditions of the massive scheme (|3^ ) to the effective 
Hamiltonia n (29|) one obtains the two- loop RG functions, that in the replica limit 
n = read 



124. 



/3„ = -(4 - d) ui 1 - i [(m+8) ii-|-12i;+2 (m+5) w] + - \2 (5m+22) u^+84:v'^+ 
6 9 L 



24 (m+b) uv+2 (14m-f 58) wu-t-4 (9to-|-33) vw+ (17to+67) 



*^ + 9 



2v^+ 
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2 (to+2) UV+ (m+2) u^ + ^^^^^w^+2 (m+2) uw+2 (m+1) vw 



«2 



(70) 



Py = -(4 -d)\v-l: [8u^+2 (m + 2) uv+2uw+2 (m + 1) vw+3w'^] + ^ 
I 6 9 

24 (m + 2) uu^+2 (3 m+6) t;u^+2 (6 to+24) uww+2 (12 m + 12) lyu 

1 2 

(3m+45) (m+8) 4u^w+(3m+9) 



44i;'' 



9 



2v*+ (m+2)wM^+ 



/ \ 2 I '^2 / \ / \ 2 

2 (m+2) H — ^ — w v+2 [m+2) uvw+2 v w 



«2 



(71) 



/3„ = -(4 - (i) 1 - i [(m + 4) w + 12v + 4w] + ^ [(5m + 27) w^+84i;^ + 
6 9 L 



4 (6m+15) vw+2 (5m+22) uw+4 (3m+18) uv+2 (m+6) 
2 (m,+2)uu+ (m+2)u^H w'^+2 [m+ 2)uw+2 (m+l)uu; 



l2 



2v' 



(72) 



70 = 



(4-rf) 



2w^ + 2 (m + 2)uv + (m + 2)u^ + 2[m + 2)uw + 



m + 3 



702 = 



2(m + l)z;w)i2, (73) 
< ^ (^2u + (m + 2) M + (m + 1) - [2v'^ + 2 (m + 2) uu+ 



(m + 2) + 2 (m + 2) uw H + 2 (m + 1) vw)ii \ , 



(74) 



Here, u,v,w are the renormalized couplings, ii,Z2 are the two-loop integrals (c. f. 
formulas (|5^), ( |53| ) and comments below them). 

Again, in a first step of the analysis one may use the e-expansion. The first order 
£-expansion results were given in Ref. iB. They can be reproduced from formulas 

@ - Pl 

in Table 



■ Tn particular one gets eight fixed points with the coordinates given 
123. T he s econd order contributions in e to the fixed poi nt c oordinates 
calculated in Ref. llHil do not change qualitatively the results obtained 1123. The main 



question of interest now is whether the above described picture of the runaway 
solution is not an artifact of an e-expansion. To shed light on this question one can 
use a more refined analysis of the fixed points and their stability, similar as it was 



done in the subsections 3.1 



6.2 



All fixed points with u > 0,v > 0,w < appear to be unstable for e > except 
of the "polymer" 0{n = 0) fixed point III which is stable for all m. However, the 
presence of a stable fixed point is not a sufficient condition for the 2nd order phase 
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Table 6. Fixed points of the RAM with isotropic distribution of the random anisotropy 
axis in e-expansion. Here, x± = (m — 2 it (m — 2)^ + 48)/8, y± = (m — 2 — 2mzi: 



^{m - 2 - 2 mzY + 4(12 - 82))/8, z = (m + 6)/(m + 8) 





u* 


V* 


W* 


I. 

II. 

III. 

IV. 

V. 

VI. 

VII. 

VIII. 









3(4-m) _ 







6 £ 


2{m-l) ^ 





S(m-l)^ 


ox _ 


12a:++m+4 

6 £ 


12a;_+m+4^ 
6m + 


12a;_ +m+4 

£ 


12H+-4z+m+4 


12H+-4z+m+4 


12H+-4z+m+4 


12j;_ ~4z+m+4 


12j;„ -4z + m+4^ 


12i/_ -4z+m+4 



transition. The fixed point should be accessible from the initial values of couplings 
and it is not the case for the location of fixed points shown in Fig. |ll]. Recall that 
the initial conditions of the RAM effective Hamiltonian in our case are given by 
uq > 0, vq > 0, Wo < and wq/vq — —m (see subsection 4.3). Indeed, starting from 
the region of initial conditions (denoted by a cross in the figure) for zero value of 
u one meets a separatrix joining the unstable fixed points I and VI and the system 
cannot reach the stable fixed point III. As far as both fixed points I and VI are 
strongly unstable with respect to u, the fixed point III is not accessible for arbitrary 
positive u either. Finally one concludes that the 2nd order phase transition is absent 
in the model as shown by runaway solutions of the RG equations. 



Ill 



X 



w 



III 



IV • 

^ VIII 



VI 



Fig. 11. Fixed points of the RAM with isotropic distribution of a local anisotropy axis. The 
only fixed points located in the octant u>0,v>0,w<0 are shown. The filled box shows the 
stable fixed point, the cross denotes typical initial values of the couplings. 



Considering the two-loop RG functions (|70|)-(|7^) directly at d = 3 and applying 
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the Pade-Borel resummation (|45|)-([47|) one gets LfJ^tiJ the fixed point coordinates 



given in the Table ^ E£3. Resummed two- loop results qualitatively confirm the pic- 
ture obtained in the first order in e-expansion: the stability of the fixedooints does 
not change after resummation. This supports a conjecture of Aharony Ea about the 
absence of an accessible stable fixed point for the RAM with isotropic distribution 
of the local anisotropy axis. 

In the other fixed points one may recover the two-loop results for the 0{m) 
model (fixed point II), the polymer 0{n = 0) model (III) and the diluted Stanley 
model (IV). The fixed point VIII contains all three couplings but is both unstable 
and not accessible from the initial values of couplings. The values of the correlation 
length and pair correlation function critical exponents v and rj, resummed in a 
similar way, are given in the table. When calculated in unstable fixed points, they 
are considered as effective exponents. 



Table 7. Resummed values of A 
two-loop approximation for d = 3 l£. 



fitted points and critical exponents for isotropic case in 



FP 


m 


u* 


V* 


w* 


V 




I 


Vm 















II 


2 


0.9107 








0.663 


0.027 




3 


0.8102 








0.693 


0.027 




4 


0.7275 








0.724 


0.027 


III 


Vm 





1.1857 





0.590 


0.023 


IV 


2 


0.9454 


-0.0322 





0.668 


0.027 




3 


0.6460 


0.1733 





0.659 


0.027 




4 


0.4851 


0.2867 





0.653 


0.028 


VI 


2 





1.4650 


-1.6278 


0.449 


-0.028 


VIII 


2 


0.7991 


0.7311 


-0.5360 


0.618 


0.026 




3 


0.8450 


0.5934 


-0.4506 


0.614 


0.025 




4 


0.8762 


0.5060 


-0.3956 


0.612 


0.025 



6.3.2. A cubic distribution of the random anisotropy axis 

Now we turn our attention to the case of the RAM with a cubic distribution of 
the random anisotropy axis (|l^). The RG functions that correspond to the effective 
Hamiltonian (|2^) in two-loop approximation were obtained in Ref. Il23 within the 
massive renormalization scheme (p2|). In the replica limit n = they read: 



Pu -£M< 1- ^ [(m+8)u+12i;+4u;+6y] +i 2 (5to+22) ^^+8 (3m+15) to+ 



8Av^+12w^+68uw+72uy+18y^+72vw+108vy+36wy 



2 (m ^- 2) uv+2v'^ +2w'^ +6uw+Av'w+6vy+6uy+6wy+3y'^ i2 



(75) 
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(3v — —el V r8u^+2 (m+2) uv+2uw+4:V'w+6vy \ + — 

6 9 



44^;■^+48^; 



12w'^v+24 (m + 2) (3to+6) vu'^+Aw^u+Au'^w+72v^y+18y'^v+ 



60uvw+36uvy+36vwy 



2u^+2 (m+2) uv^+ (m+2) vu'^+ 



2w v+6uvw+Av w+6v y+6uvy+6vwy+3y v 12 



(76) 



= -ewJ 1 - 1 [8w+12z;+4u+6y] + ^ 44w2+84t)2+120'u;w+2 (m+6) ^2+ 



68m«7+72w?/+182/2+2 (6m + 36) ™+108t;y+36-uj/ 



n- 



2 r 
9 



(m+2)u^ + 



2 (m + 2) uv+2v^ +2w'^ +&uwA-Avw+&vy+&uy+&'wy+'iy'^ 



«2 



(77) 



(3y = -e\ y- \ [9y^+8ww+12wy+12uy+12wy] +\ (4m + 72) u^w+72w^u+ 
16 9 L 

54y^+84i;^y+ (6m + 84) u^j/+84w^?;+144j/^i;+144y^u+144y^u;+ 



96tiww+2 (6m + 84) uvy+2b2uwy+lQ8vwy 



*^+9 



2w^?;+2 (m + 2)M'yy+ 



(m+2) u y+2w y+Guwy+Avwy+dy v+6y u+6y w+3y 



«2 



(78) 



70 = -- {2v^ + 2 (m + 2)ut; + (m + 2)^^ + 6uw + 2^^ + Avw + 6wy+ 
9 

6uj/ + 6w7y + 37/^) 12, (79) 
702 = I i i ^2v + {m + 2)u + 2w + y^ - (2w^ + 2 (m + 2) + (m + 2) m^+ 



6uii; + 2w^ + Avw + 6u?/ + 6uy + 6wy + 3?;^) ii > , 



(80) 



Again, ii, 12 stand for the two-loop integrals and the renormalized couplings u, v, 
w and y correspond to the effective Hamiltonian ( p9| ) . 

As in the isotropic case the first RG result for the RAM with cubic distribution 
of local anisotropy axis was obtained in the first order in e c3. It brought about 
the evidence of 14 fixed points and may be easily reproduced from the functions 
(|7^)-(|80|) removing the two-loop contributions. The coordinates of the fixed points 
are listed in the Table | @ 

The results of the linear e analysis El states that among the fixed points with 
u > 0,v > OjW < only a "polymer" 0{n = 0) fixed point III is stable for all m for 
£ > 0, but it is not reachable from the initial values of couplings (see Fig. |l^). The 
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Table 8. 
Here, q± 



Fixed points of the RAM with the random cubic anisotropy distribution (e-expansion) . 
= (m - 



4± Vm2 + 48)/8, /3± = -(m + 12± + 48)/6, A±± = 6o± + 3/3± + m + 6. 

The fixed points XV-XVII appear only in the two-loop approximation due to the degeneracy of 
the corresponding one-loop functions. 



I. 

II. 

III. 

IV. 

V. 

VI. 

VII. 

VIII. 

IX. m ^ 2 

X. 

XI. 

XII. 

XIII. 

XIV. 



XV. 
XVI. 

XVII.m = 2 









6 

4(m-l) 












6(m-4) 
16{m-l) ' 

1 = 



m — 4 ^ 
4(m-2) ^ 



3a I 

— ^ 

3a „ 









3(m-|-4) 



4Aj 



+ 



3{m-|-4) 



4A 



3(m-|-4) 
4A_+ ^ 
3(m-|-4) 

iA__ ^ 








y 








9^ 




2(m-4) 

3m, ' 
m — 4 , 
3(m-2) ' 

3(3+ 



3/3+ 

-4- + 

3/3- 




reason is a separatrix joining the unstable fixed points I and VII and separating 
initial values of couplings (shown by a cross in the Fig. |l^ and fixed point III. 
The possible runaway behavior of the RG flow lead tp_the conclusion about a 
smearing out of the phase transitinn as Tc is approached 112. However, a subsequent 



study of Mukaniel and Grinstein li^ brought about the possibility of a second order 
phase transition with the scenario of a weakly diluted quenched Ising model 119. 
Indeed, accounting for the terms qualitatixehLxhanges the picture: performing 
the perturbation theory expansion one gets iHl-llHa not only the corrections to the 
coordinates of the fixed points I-XIV (we do not display them here as far as they as 
are too cumbersome) but also the new fixed points XV, XVI, XVII (see the bottom 
of the Table ||). The appearance of the pairs of fixed points denoted by XV and 
XVI is caused by the fact that the /3-functions Pw, f3y at u = v = (/?„, f3y at 
w — y = 0, correspondingly) are degenerated at the one loop level. This leads to 



the •\/e--expansion (c.f. subsection 6.1) of the weakly diluted quenched Ising model. 
The ^/e expansion of the fixed point XVII holds only for to = 2 and is caused by 
the one-loop degeneracy of the Pv, Py functkms for w = (c. f. singularity at 
TO 2 in the e-expansion of the fixed point IX) E£3. 

Checking the stability of the new fixed points one finds that all of them are un- 
stable except the Ising-likc fixed point with w < 0, y > of the pair XV. Moreover, 
this point is reachable from the initial values of the couplings. As far as this is the 
fixed point of the diluted Ising model one concludes, that in the critical region the 
RAM with cubic distribution of the random anisotropy axis ( |l5| ) decouples into to 
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independent dilute Ising models and the phase transition for any m is governed by 
the familiar random Ising model critical exponents. However, it is worthwhile to 
keep in mind that above picture is obtained in the frames of the "naive" analysis of 
the e (and -i/e) expansion and it is highly desirable to confirm it by a more reliable 
analysis of fixed points and their stability. 



w 




Fig. 12. Fixed points of the RAM with distribution of a local anisotropy axis along hypercube 
axis for u = 0. The only fixed points located in the region v > OyW < are shown. Filled boxes 
show the stable fixed points, the cross denotes typical initial values of couplings. 



The results of the application of the Pade-Borel resummation nrn ced ure (Ba) 



) to the two-loop /3-functions (|7q)-(|78|) at d = 3 are given in Refs. LfJ'EJfa where 16 



fixed points were obtained for this model. In the Table |^ we present the numerical 
values of fixed points coordinates with u* > 0, v* > 0, w* < 0. The last fixed 
point XV in Table ^ corresponds to the stable fixed point of y^-expansion. It 
has coordinates with u* = v* = 0, w* < and y* > and is accessible from the 
typical initial values of couplings (marked by a cross in the Fig. 12). The Pade-Borel 
resummation procedure (|75|)-(|78|) does not reveal any other stable fixed points in 
the region of interest. Thus one draws the conclusion that the system in the critical 
regime is characterized by the same set of critical exponents as the weakly diluted 
quenched Ising model. 

In the other fixed points, one recovers the familiar two- loop numerical results 
for the Gaussian (fixed points I, VII), m-vector (II), polymer 0{n = 0) (III), Ising 
(V, X), diluted m-vector (VI), and cubic (VIII) models. Fixed point IX belongs to 
a new universality class. In the Table ^ we give the numerical values of the critical 
exponents in these fixed points as well: if the flow from the initial values of couplings 
passes near these fixed points one may observe an effective critical behavior governed 
by these critical exponents. Note that the fixed point VI has a singularity for to = 1 
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which formally leads to the ^/e expansion. However, as it follows already from the 
microscopic Hamiltonian (^3|) for m ~ 1 the RAM degenerates to the regular Ising 
model. This also may be seen from the effective Hamiltonians (|2^), (p9|): in the 
replica limit at m = 1 v and w terms (and for ( p9| ) also u and y terms) appear to 
be of the same symmetry. Furthermore, couplings v and w appear to have opposite 
signs and the same absolute value, therefore they are can celed and the problem 
reduces to the scalar 0"* theory with only one coupling t22 . 



As a possible generalization of the RAM one may consider the case when the 
quenched randomness is present in both the random-site and the random-anisotropy 



form. This is described by effective Hamiltonians (p8D, (29) where the coupling ?;o 
may be of either sign. We have checked the region w < for the presence of new 
fixed points and verified that they are absent. Therefore, again the fixed point XV 
is the only reachable and stable one and the observed critical behavior is unique. 
As an interesting extension of these studies one may consider the influence of other 
types of disorder in combination with the random anisotropy as well as of other 
random axis distributions. 



Table 9. Resummed values of tbe-fi*6d points and critical exponents for cubic distribution in 
two-loop approximation for d = 3 t23'E£f|. 



FP 


m 


u* 


V* 


w* 


y* 


V 




I 


Vm 














1/2 







2 


0.9107 











0.663 


0.027 


II 


3 


0.8102 











0.693 


0.027 




4 


0.7275 











0.720 


0.026 


III 


Vm 





1.1857 








0.590 


0.023 


V 


Vm 











1.0339 


0.628 


0.026 


VI 


3 


0.6460 


0.1733 








0.659 


0.027 




4 


0.4851 


0.2867 








0.653 


0.027 


VII 


Vm 





2.1112 


-2.1112 





1/2 







2 


1.5508 








-1.0339 


0.628 


0.026 


VIII 


3 


0.8393 








-0.0485 


0.693 


0.027 




4 


0.5259 








0.3624 


0.709 


0.026 


IX 


3 


0.7096 


0.1695 





-0.1022 


0.659 


0.027 




4 


0.4190 


0.2751 





0.1432 


0.653 


0.027 


X 


Vm 





0.6678 


-0.6678 


1.0339 


0.628 


0.026 


XV 


Vm 








-0.4401 


1.5933 


0.676 


0.031 



As we have seen from the analysis of the subsection 6.3 the local anisotropy axis 
distribution has a crucial influence on the scenarios of low-temperature behavior. 
For the RAM with isotropic distribution the inaccessibility of any stable fixed point 
leads to the conclusion that a second order phase transition is absent. The existence 
of a reachable disordered Ising-like fixed point in case of the RAM with cubic distri- 
bution of the random anisotropy axis indicates that the system undergoes a second 
order phase transition with the critical exponents of the weakly diluted quenched 
Ising model. Here, we want to attract thereaders attention to the similarity of the 
critical behavior of both the random-site IIj and the random-anisotropy lifl quenched 
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magnets: if new critical behavior is present then it is always governed by the critical 
exponents of the site-diluted Ising type. The above calculations of a "phase dia- 
gram" of the RAM are based on two-loop expansions improved by a resummation 
technique. Once the qualitative picture becomes clear there is no need to go into 
higher orders of a perturbation theory as far as the critical exponents of site-diluted 
Ising model are known by now with high accuracy (see Table |^ and discussion of 



the subsection 6.1) 



7. Conclusions 

To conclude this review, let us return to the question posed in the introduction: 
what information about the influence of weak quenched disorder on criticality can 
one get analyzing the RG expansions? Comparing the results of the theoretical 
analysis in Section ^ with the available data about different physical phenomena 
listed in Section ^ one definitely can state that in most cases the RG analysis 
provides a correct description. Indeed, the change of the universality class of the 
Ising model upon weak quenched dilution predicted by the RG analysis is confirmed 
by numerous experiments on diluted uniaxial magnets FepZni_pF2, MnpZni_pF2 as 
well as by extensive MC simulations (c. f. data of the Table || and of the Section ||). 
By high precision experiments under microgravity conditions it was shown c3 that 
the specific heat asymptotic critical exponent of liquid He^ for the A-transition is not 
influenced by quenched disorder and belongs to the universality class of the m = 2 
Stanley model. This result also follows from the RG analysis as far as the Stanley 



131 



model heat capacity critical exponent a changes its sign already fnr t" < 2 
The MC simulations of self-ampiding walks on random-site lattices EaE2 confirmed 
the theoretical result of Ref. e3 that universal properties of self-avoiding walks are 
not influenced by point-like uncorrelated disorder. The change of the universality 
class of self^voiding walks by a long-range-correlated disorder predicted by the RG 
analysis EJ'Ej waits for an experimental or simulational confirmation. The MC study 
of the influence of the long-range-correlated disorder on the ferromagnetic phase 
transition confirms the change of the universality class in this case E3. One more 
interesting check of the RG predictions might be the verification (experimental 
or in MC simulations) of the envisaged random Ising model critical behavior for 



the 



127 



128 



12£ 



anisotropy magnet with a cubic distribution of the anisotropy axis 

132f 



As we have seen in these lectures, for the majority of the phenomena listed above 
a careful analysis of the RG expansions provides also a quantitative description, 
giving reliable numbers for the physical values describing criticality. However, if 
one compares the RG expansions for non-diluted ( "pure" ) and diluted systems one 
certainly notices essential differences. A typical scenario that one encounters when 
analyzing the critical properties of the d = 3 Stanley model (i.e. applying the field 
theoretical RG technique to the 0(TO)-symmetrical (j)^ model) may be sketched by 
two steps: (i) already in the first non-trivial order of the perturbation theory one 
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observes a new non-trivial critical behavior (a stable non-trivial fixed point); (ii) 
the higher order contributions do not change the picture qualitatively and improve 
the numerical outcome provided the resummation is applied. This scenario does not 
hold for most of the models with disorder. Here, the first-order perturbation theory 
expansions as a rule do not lead to new criticality (a non-trivial fixed point is absent 
or its stability is questioned). The (experimentally confirmed!) answer about new 
critical behavior is obtained only by considering higher order contributions with 
the help of the resummation technique. But then again, even the resummed series 
appear to be divergent starting from some order of perturbation theory showing an 
"optimal truncation" behavior. Some expansions that are familiar and useful for 
the "pure" models appear to be of no use for certain disordered models: e.g. the 
e-expansion turns into the -y/e-expansion for the random-site Ising model and does 
not lead to any quantitative information. 

Technically, this "strange" behavior of the perturbation theory expansions for 
the disordered models is not strange: it rather signals the diverg ent nature of the 
series. For some disordered models this behavior was shown to originate from 
Griffith-like singularities H caused by zeroes of the partition function of the pure 
system. The physical reason for Griffith singularities is the existence of ordered 
islands within the temperature interval below the critical temperature of the pure 
model and above that of the diluted model. So in this sense the studies of random 
models reviewed above may serve as examples when more complicated physics leads 
to a more complicated theoretical description. And as one may conclude from the 
results obtained so far in many cases this description is a fruitful one. 
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